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FOREWORD 


This  report  summarizes  the  results  of  a study  performed  by  Caywood-Schilkr  Division 
of  A.T.  Kearney,  Inc.  under  U.S.  Air  I;orce  Contract  F33615-73-C-2Q78.  The  work  was 
conducted  between  1 July  1973  and  31  March  1974,  under  the  direction  of  the  Air  Force 
Aero  Propulsion  Laboratory,  with  Mr.  G.  W.  Gander  (AFAPL/SFH)  acting  as  Project 
Engineer. 

Work  was  sponsored  by  JTCG/AS  as  part  of  the  3-year  TEAS  (Test  and  Evaluation, 
Aircraft  Survivability)  program.  The  TEAS  program  was  funded  by  DDR&E/ODDT^IL  The 
effort  was  conducted  under  the  direction  of  tire  JTCG/AS  Vulnerability  Assessment  Sub- 
group, as  part  of  TEAS  element  5. 1.6.6,  Development  oj  Models  for  Assessment  of  the 
Vulnerability  of  Aircraft  Fuel  Systems. 

A study  was  conducted  to  develop  a dynamic  model  of  the  vulnerability  of  an  aircraft 
fuel  system  to  threats  posed  by  hostile  weapons.  Improvement  was  achieved  in  treating  fuel 
system  vulnerability.  Further  development  of  tiie  fuel  system  model  is  recommended. 


DISCLAIMER 

Estimates  in  this  report  are  not  to  be  construed  as  an  official 
position  of  any  of  the  Services  or  of  the  Joint  AMC/NMC/ 
AFLC/AFSC  Commanders. 


NOTE 

Information  and  data  contained  in  tltis  document  are  based 
on  reports  available  at  the  time  of  preparation,  and  the  re- 
sults may  be  subject  to  change. 
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— Tliis  report  presents  the  results  of  a study  to  develop  a 
dynamic  ‘model  of  the  vulnerability  of  an  aircraft  fuel  system  to 
tiie  threats  posed  by  hostile  weapons.  A Monte  Carlo  model  was 
developed  to  calculate  the  probability  of  hit  along  segments  of  a 
specified  flight  profile  at  each  point  where  a specified  weapon 
system  could  pose  a threat  to  the  fuel  system.  An  Air  Force  devel- 
oped computer  model  (Well-Stirred  Fuel  Tank  Model)  is  used  to 
compute  fuel  state  in  each  fuel  tank  under  study  at  increments 
along  the  flight  path.  These  are  used  as  inputs  to  the  Monte  Carlo 
model. 

Given  that  a hit  takes  place,  the  probable  trajectories  (liquid- 
air,  liquid-liquid,  and  air-air)  are  calculated,  and  the  probabilities 
of  lethal  outcomes  (explosion,  internal  fire,  external  fire,  leak)  are 
computed.  The  model  ranks  the  most  likely  events,  and  a hazard 
index  is  generated  which  portrays  the  most  important  threats  to 
the  fuel  system  on  the  specified  flight  path. 

Tlie  resulting  model  gives  an  improved  measure  of  the  impart 
of  fuel  state  on  the  vulnerability  of  a fuel  system  on  aircraft 
vulnerability.  It  does  not  incorporate  consideration  of  the  effects 
of  fuel  slosh,  vibration,  or  vent  geometry.  Further  refinement  and 
development  is  recommended.  ^ 
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INTRODUCTION 

PURPOSE 

This  report  presents  tire  results  of  an  exploratory  study  to  develop  a dynamic  represen- 
tation of  tiic  vulnerability  of  an  aircraft  fuel  system  to  hostile  threats. 

BACKGROUND 

Some  analytical  studies  use  gross  aggregation  in  representing  the  vulnerability  of  an 
aircraft  fuel  system.  Typically,  a fixed  percentage  of  fuel  is  assumed,  empty  and  external 
fuel  tanks  ignored,  and  the  wide  variation  in  probability  of  a reaction  within  the  various 
kinds  of  fuel  tanks  is  compressed  to  a single  probability  of  fire  given  a hit.  Also,  distinction 
in  hazard  between  liquid  fuel  and  tank  ullages,  influence  of  tank  wall  and  liquid  temper* 
atures  on  the  ullage  composition,  tank  overpressure  caused  by  internal  fires,  and  effect  of 
venting  are  largely  ignored.  Analyses  conducted  under  these  simplifying  assumptions  arc 
valid,  and  study  results  arc  reliable  and  important;  but  more  precise  representation  is 
required  for  test  planning,  fuel  system  design,  and  detailed  examination  of  the  fuel  system 
vulnerability  area.  The  greater  accuracy  of  a dynamic  fuel  system  vulnerability  model  will  be 
beneficial  in  future  aircraft  development  studies  because  delineation  of  specific  vulnerabil- 
ities and  more  precise  measurement  of  previously  identified  vulnerability  will  be  possible 

In  performing  this  background  study,  the  WSFT  (Well-Stirred  Fuel  Tank)  computer 
program,  developed  for  the  Air  Force  by  Dynamic  Science,  Inc.,  was  considered  the  most 
accurate  representation  of  internal  fuel  state.  The  WSFT  program  docs  not  consider  the 
hostile  threats;  therefore,  a vulnerability  model  was  constructed  to  combine  the  WSFT 
program  fuel  states  with  anticipated  threats.  . 

WSFT  PROGRAM 


GENERAL 

The  WSFT  program,  which  is  an  integral  part  of  the  fuel  vulnerability  model,  was 
developed  under  a prior  AFAPL  study1 . The  program  determines  the  fuel-to-air  ratio  in  the 
ullage  space  of  a fuel  tank  as  a function  of  time  for  a particular  input  mission  profile  and 
describes  the  state  of  fuel  and  vapor  space  in  fuel  tanks,  accounting  for  mass  and  energy 
transport  due  to: 

( 1 ) fuel  evaporation 

(2)  venting  effects 

(3)  heat  transfer  between  ullage,  tank  walls,  and  liquids 

(4;  outgassing  as  dissolved  air  is  removed  from  the  liquid  as  the  aircraft  climbs. 

The  program  does  not  consider  interaction  of  a fuel  tank  and  an  ignition  source,  such  as  in 
an  incendiary  projectile.  The  fuel  vulnerability  model  developed  in  this  study  combined  the 
interaction  of  threat  and  fuel  tank. 

* Air  Force  Aero  Propulsion  Laboratory.  Analysis  of  Aircraft  Fuel  Tank  Fire  and  Explosion  Hazard , by  T.C.  Kosvic. 
L.B.  Zung,  M.  Go r stein.  Wright-Pattcrson  Air  Force  Base,  OH,  AP'APL,  March  1971.  75  pp.  (AFAPL-TR-71-7,  publication 
UNCLASSIFIED.) 
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INPUTS  RLQUIRED 

Basic  inputs  required  by  the  WSFT  program  are: 

1 . Altitude  profile  of  aircraft 

2.  Liquid  temperature  history 

3.  Skin  and  structure  temperature 

4.  Vapor  pressure  relations  of  liquid  fuels 

5.  Ullage  volume  and  exposed  surface  area  schedule 

6.  Vent  size 

7.  Internal  heat  transfer  coefficients  between  ullage  and  tank  structure. 

Standard  FORTRAN  IV  NAMELIST  is  used  for  all  input  data  to  the  WSFT  program. 
The  variable  list  and  definitions  for  the  input  data  under  the  NAMELIST  name  DATA  are 
presented  in  Table  1.  An  example  of  a set  of  inputs  to  the  program  is  provided  in  Figure  1 . 


OUTPUT 

The  normal  output  of  the  WSFT  program  is  a report  for  each  print  time  throughout 
the  mission  profile.  Each  report  describes  the  state  of  the  fuel  and  vapor  space  at  that 
particular  time  in  terms  of  air  and  fuel  partial  pressures,  fuel  vapor  pressure,  fuel-to-air  ratio, 
and  mass  and  mass  flux  of  fuel  vapor  that  has  been  vented,  evaporated,  outgassed,  or 
condensed.  In  addition  to  these  parameters,  the  altitude,  speed,  and  amount  of  fuel  used  are 
printed  out.  An  example  of  the  printout  of  this  report  is  shown  in  Figure  2.  The  WSFT 
output  parameters  used  for  the  vulnerability  model  were  the  fuel-to-air  ratio  and  the 
amount  of  fuel  used  as  a function  of  time.  Tables  of  fuel-to-air  ratios  and  the  percent  of  fuel 
remaining  in  each  tank  as  a function  of  time  were  generated.  These  tables  were  stored  on 
magnetic  tape  and  used  as  input  to  the  fuel  vulnerability  model. 


INHERENT  ASSUMPTIONS  AND  LIMITATIONS 

The  WSFT  program  assumes  a homogenous  mixture  of  fuel  vapor  and  air.  The  mixing 
of  air  and  vapor  is  assumed  to  occur  rapidly  with  no  appreciable  difference  in  fuel-to-air 
ratio  within  the  ullage  volume.  The  program  is  particularly  applicable  for  shallow  tanks  or 
tanks  where  he  ratio  of  ullage  volume  to  liquid  fuel  surface  is  small.  It  was  considered 
beyond  the  scope  of  this  study  to  develop  a new  model  which  would  incorporate  the 
concept  of  a stratified  ullage  space  of  different  fuel-to-air  ratios. 

The  effects  of  vibration  and  slosh  on  the  fuel-to-air  ratio  are  not  included  in  the  WSFT 
program.  One  of  the  inputs  required  is  the  liquid  temperature  history  of  the  fuel  in  the  tank. 
This  information  must  be  provided  by  the  user;  it  is  important  because  aerodynamic  heating 
or  cooling  may  have  a significant  influence  on  the  temperature  of  the  fuel  and  heat  sources 
within  the  aircraft.  These  factors  should  be  considered  in  the  development  of  more  sophis- 
ticated models  for  determining  fuel-to-air  ratios. 
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Table  1.  Variable  Names  and  Definitions. 


Name 

Definition 

Units 

SDATA 

Ident  of  data  block 

RGAS 

Universal  gas  constant 

ft-lbm/lb-mole  °R 

EMWA 

Mass  of  air 

ibm/lbm-mole 

EMWF 

Mass  of  fuel 

lbm/lbm-molc 

CPA 

Specific  heat  of  air 

BTU/lbm  °R 

CPF 

Specific  heat  of  fuel 

BTU/lbm  °R 

TA 

Temperature  of  the  ullage 

°F 

HJ(J) 

Hj  heat  transfer  film 
coefficient 

BTU/ft2  hr°R 

J=l,  to  fuel  surface 

J=2,  to  side  of  tank 

J=3,  to  top  of  tank 

ZOG 

1 Set  equal  to  one 

None 

DV 

Diffusion  coefficient 

ft2/hr 

CDELTA 

Characteristic  length  for 
evaporation 

ft 

KTANK 

0 Set  equal  to  zero 

None 

GALO 

Initial  volume  of  fuel 

gal 

BETA 

Bunsen  coefficient 

CONI 

Outgassing  coefficient 

hr*l 

CON  2 

Solution  coefficient 

hr1 

TVENT 

If  TVENT=0,  incoming  air  will  be 
calculated  from  altitude 
and  Mach  number  schedule 

°F 

If  TVENT^O,  all  incoming  air 
will  have  temperature=TVENT 

TF(1,1) 

' 

Table  of  time  values 
corresponding  to  fuel 
temperature  table 

hrs 

/ 
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Table  I.  Variable  Names  ami  Definitions  (t  ontinued). 


Name 

Definition  ! 

Units 

TNI  ,2) 

Table  of  fuel  temperatures  at 
times  corresponding  to  TF(  t ,1 ) 

"!■ 

TSIDE(  1,1) 

Table  of  time  values 
corresponding  to  tank  side 
temperature  table 

hrs 

TS1DL(1,2) 

Table  of  tank  side  temperatures 
at  times  corresponding  to  TSIl)h(  1,1) 

°F 

TTOP(l,l) 

TTOP(l,2) 

Same  as  TSiDF  except  applies 
to  top  of  the  tank 

ALT(l.l) 

Table  of  times  corresponding 
to  altitude  schedule  tabic 

hrs 

ALTO  ,2)  ' 

Table  of  altitudes  at  times 
corresponding  to  ALT(ltl) 

KFT 

GALDOT(l,l) 

Table  of  timp  values 
corresi  Hiding  to  fuel  usage 

hrs 

GALDOT(l,2) 

Table  of  fuel  usage  at  times 
corresponding  to  GALDOTfl  ,1 ) 

gal/hr 

EMINF(I.l) 

Table  of  time  values 
corresponding  to  flight 
Mach  number  schedule 

hrs 

EMINF(1 ,2) 

Table  of  flight  Mach  numbers 
at  times  corresponding 
to  FMINF(l.l) 

None 

PVAP(1,I) 

Table  of  temperatures  corresponding 
to  fuel  vapor  pressure 

Op 

PVAP(1,2) 

Table  of  fuel  vapor  pressure 
corresponding  to  PVAP(1,1) 

psia 

TO 

initial  time  for  start  of 
integration 

hrs 

TMAX 

Final  integration  time 

hrs 

DT 

Time  step  for  integration 

hrs 

DTPRNT 

Print  time  interval 

hrs 

AV 

Area  of  the  vent 

ft  2 

DELHF 

Heat  of  formation  of  fuel 

BTU/lbm 
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1 able  1.  Variable  Names  and  Definitions  (Continued). 


Name 

Definition 

Units 

— < — — " — 

ULLGH 

Ullage  length 

ft 

ULW1D 

Ullage  width 

ft 

ULMT 

Ullage  height 

ft 

S 

End 

SDATA 

RGAS=1545, 

EMWA=28.966,EMWF=72, 

CPA=0.24,CPF=0.49, 

TA-60, 

HJ=3*2 

U L LG  H= 1 0.0,  U LW I D= 1 0.0,  U L HT=0.55, 

DELHF=1, 

ZOG=1, 

AV=0.16, 

DV-0.3, 

CDELTA=0.01, 

KTANK=0, 

GALO=5573. 

BETA=0.16, 

CO  N 1 = 1 OOO.CO  N2=0, 

TVENT=70, 

TF<1.1)=0,1, 2, 3,4.5,6,7,8,9,10,1 1,12, 

TF  ( 1 ,2)=60,60, 58.48.48, 45,1 5,1 8,20,25,1 10,1 30,1 20, 

TSIDEO.I  )=0,12,TSIDE(1,2)=70,70, 

TTOP(  1 , 1 )=0, 1 2,TTOP(  1 ,2)=70,70, 

ALT(1 ,1  )=0,.1  ,.3,.5,1 ,4,5,6, 8, 8.3,1 2, 

ALT(  1 ,2)=0, 8,20,22,22,20,1 8, 20, 22, .25, .25, 

GALDOTO  .1 ) =0,8.3, 8.31 ,10.3,10.31 ,1 2, 

GALDOTd  ,2)=0,0,2779, 2779,0,0, 

PV  AP(1,1)=17,41,67 ,96,129,166, PVAP(1,2)=.35,.60,1.1. 2.0, 4.0.8.0, 

EMINF(1,1)=0,.5.12  ,EMINF<1,2)=0,.85,.85, 

TO=0,TM  AX=1 1 ..DTP 8 NT"  .2r;,DT- 001$ 


*A  header  card  must  be  present  for  each  run  on  the  WSFT  program.  The  subsequent  data 
cards  contain  the  variable  names  and  corresponding  values  required  by  the  model.  The  data 
is  punched  on  cards  beginning  in  column  two. 

Figure  1 . Example  of  a Set  of  Inputs  for  WSFT  Program. 
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Time  hrs  = 1.7500  Air  partial  pressure  = .77882E+03 

Mach  number  =.85000  Fuel  partial  pressure  = ,13494E+03 

Vent  velocity  = 1.6085  (ft/hr)(positive  into  tank)  Fuel  vapor  pi  essure  = .13486 E+03 

Integration  error  =-1.6841  Air-fuel  ratio  = .23219E+01 

total  mass-percent 
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VULNERABILITY  MODEL 


LTllLli  Y 

Ylu:  vulnerability  model  was  designed  to  ascertain  the  most  hazardous  phases  of  a 
designated  flight  protile,  given  a specific  aircratt  and  a specified  set  of  hostile  weapons.  It  is 
intended  tor  use  in  support  of  laboratory  'esfing,  both  to  reduce  the  number  of  tests 
required  and  to  eliminate  unnecessary  tests. 


INPUT'S 

lo  prepare  for  model  runs,  it  is  necessary  to  iiave  input  information  ;n  several  catego- 
ries (i-e..  mission,  weapons,  and  aircraft).  See  the  Appendix  for  program  fisting  and  sample 
case. 

In  addition  to  the  input  from  the  WSR  program,  the  vulnerability  model  requires 
input  information  regarding  ignition  and  detonation  probabilities  as  functions  of  residual 
IK’ne  tutor  energy,  and  hydraulic  ram  probability  as  a function  of  impact  kinetic  energy.  It 
also  requires  a geometric  description  ot  the  fuel  system  and  the  hostile  weapons  under 
consideration. 

■Miaittn 

The  flight  profile  must  be  selected;  which  includes  altitude,  velocity,  and  time  of  the 
target  aircraft  during  the  period  when  Hostile  weapons  may  be  expected. 

Weapons 

The  model  wili  accept  any  mixture  of  kinetic  energy  pc  net  tutors  or  fragmenting  war- 
h."  d weapons  in  a sector  of  attack  determined  by  a range  of  permissible  azimuth  and 
elevation  cnglcs  for  each  hostile  weapon  system. 

KINETIC  ENERGY  WEAPONS.  The  parameters  • diich  characterize  the  kinetic  energy 
penetrators  (e.g..  ball,  AP*.  and  API  ) are;  mass  of  the  penetrator,  time  of  incendiary 
ignition  and  incendiary  burnout  (relative  to  initial  contact),  muzzle  velocity,  drag 
coefficient,  and  a table  of  mil  aiming  errors  as  a function  of  target  velocity.  For  each  shot, 
the  azimuth  and  elevation  angles  are  taken  from  uniform  distributions  within  the  prescribed 
limits  for  the  firing  weapon  system.  The  probable  aiming  error  is  calculated  from  the  table 
of  mil  errors.  A particular  Dm  (miss  distance)  is  chosen  from  a normal  distribution 
characterized  by  this  probable  error.  A point  is  chosen  at  random  on  the  circumference  of  a 
circle  having  a radius  equal  to  the  Dm  and  lying  in  the  plane  perpendicular  to  the  relative 
velocity  vector.  This  point  and  the  relative  velocity  vector  determine  the  trajectory.  If  the 
trajectory  intersects  any  of  the  fuel  tanks,  a hit  is  said  to  occur  on  that  tank.  If  a hit  takes 
place,  the  geometry  of  the  situation  and  the  conditions  within  the  fuel  tank  determine 
whether  a particular  damage  mechanism  occurs. 


Armor-picrcing  ami  armor-piercing  incendiary. 
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FRAGMENTING  WARHEAD  WEAPONS,  The  parameters  which  characterize  the  frag- 
menting warhead  weapons  are:  (1>  average  fragment  mass,  < 2)  fragment  static  emission 
, velocity.  U)  fragment  slowdown  constant,  (4)  static  fragment  spray  limits,  (S)  average  static  <. 
fragment  dcority/^teradiac,  (6)  rauzzk  velocity  of  projectile,  (7)  Slowdown  constant  of  pro-  v 
Jectik,  f S)  aimmg  rigma,  and  fazing  sigma.  For  each  shot,  a trajectory  is  chosen  as  with 

theisn&d  shot  weapons,  using  the  aiming  sigma  in  place  of  the  table  of  mil  errors;  A burst' 
point  is  chosen  along  the  trajectory  from  a normal  distribution  using  the  fuzing  sigma.  For 
each  foci  tank,  the  riaveRhpe  equation  of  the  fragments  iv  solved  using  an  iterative  proce- 
V dure:  The  expected  number  of  hits  on  each  tank  is  calculated;  this  number  and  the  coodi- 
. tioosin  the  tank  determine,  the  probability  of  a particular  damage  mechanism’  '■"}/: ’■  yb  :: 

I /.'V.  ■ ' : r>.’  " 

Aircraft'  >'•'<>"  • ;• 

rrrfrr,  vv.*y;.v  ■••••. -•  <b-b  ; ( ...  y.t';7:y;.  y^p.i  u 

' The  tarset  aircraft  is  represented  as  a set  of  fuel  tanks.  Each  tank,  is  a rectangular 
parallelepiped  characterized  by  *he  coordinates  of  its  center,  and  by  its  length,  width,  and 
bright  Constats,  which  must  be  Supplied  by  the  user,  are  stored;  they  represent,  for  the  : 
aspect  of  each  tank,  whether  an  external  ignition  source  is  present,  and  the  distance  of  that 
aspect  from  the  skin  of  the  aircraft.  The  target  is  further  characterized  by  the  target  velocity 
and  by  the  energy  levels  required  to  produce  penetration  of  each  fuel  tank  and  hydraulic 

;y;.  '.fb  'V'"^ 

t by*-,  ■'['  r y ' ' '’’V  -v^  ■ 

PROGRAM  OUTPUT  bJby  ' • s ' ’..y"  ^ 

There  are  t wo  reports  generated  by  the  program.  Report  1 is  a summary  of  fuel  states 
for  each  tank  as  a function  of  time.  (See  Table  2).  Report  2 is  a summary  of  each  non-zero 
hazard  incident.  (See  Table  3.)  It  snows  the  fuel  state  at  the  time  of  the  incident,  summa- 
, nzes  the  probable  trajectories  through  the  tank  if  impact  occurs,  and  presents  probabilities 
of  no  effect,  leak,  external  fire,  instructive  ram,  and  internal  fire /ex  plosion.  These  proba-  ... 
briities  are  calculated  as  the  average  over  several  Monte  Carlo  trials  for  each  combination  of 
fuel  tank,  weapon,  and  time  intervals  resulting  in  a non-zero  Ph  (probability  of  a hit)  oh  the  : 
fuel  tank.  A summary  of  all  Report  2 outputs  is  made,  arranging  hazards  in  descending 
'order,  ■■■..'  ' ''  .7  \ v ",  7 ’ . r V: ; 

For  solid  shot  weapons,  the  probabilities  of  liquid  and  vapor  exit  are  given.  On  the1 
basis  of  ail  these  probabilities,  a hazard  index  is  calculated;  which  indicates,  on  a scale  of 
zero  through  one,  the  relative  likelihood  of  lethal  damage  occurring  to  the  aircraft  as  a result 
of  this  encounter.  .. 

COORBlNATt  SYSTEM  AN DAZ1MUTHE OVATION  CONVENTION 

The  coordinate  system  used  in  this  model  has  its  origin  at  the  cen  ter  of  gravity  of  the  . 
target  aircraft;  therefore,  it  is  a moving  coordinate  system.  The  X axis  is  positive  in  the 
direction  of  travel.  (See  Figure  3.)  The  Y axis  is  positive  in  the  direction  of  the  left  wing, 
and  the  2 axis  is  positive  in  an  upward  direction. 

In  the  analysis,  the  orientation  of  certain  vectors  with  respect  to  certain  axes  is  some- 
times  expressed  in  terms  of  direction  cosines  and  in  terms  of  azimuth-elevation;  thus, 
conversion  between  these  terms  is  required.  The  sign  convention  for  azimuth-elevation  needs 
definition  because  all  users  do  not  use  identical  conventions. 
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Figure  3.  Moving  Coordinate  System. 


Consider  the  X,  Y.  and  Z axes  and  the  Vsh  (shell  velocity)  having  direction  cosines 
^X’  My>  an(J  l*Z'  The  ^ (elevation  angle)  is  the  angle  between  the  vector  and  the  XY  plane. 
It  is  considered  positive  if  the  vector  has  an  upward  component,  and  negative  if  it  has  a 
downward  component.  The  Az  (azimuth  angle)  is  the  angle  between  the  projection  of  the 
vector  on  the  XY  plane  and  the  negative  X axis.  It  is  considered  positive  if  the  vector  has  a 
component  in  the  -Y  direction. 

The  relations^between  the  direction  cosines  and  azimuth-elevation  can  be  derived  by 
resolving  the  unit  VgH  into  its  components  along  the  coordinate  axes.  The  magnitude  of  the 
Z axis  component  is  |sin  El  |,  while  the  projection  in  the  XY  plane  is  |cos  El  J.  The  latter 
component  may  be  projected  on  the  X axis  to  yield  cos  El  cos  Az  and  on  the  Y axis  to 
yield  I cos  El  sin  Az  |.  Thus,  the  unit  vector  has  components  along  the  X,  Y,  and  Z axes 
whose  magnitudes  are  |cos  El  cos  Az  j,  |cos  El  sin  Az|,  and  | sin  El  |.  Taking  the  sign 
conventions  into  account: 

= -cos  El  cos  Az  (1) 

Py  = -cos  El  sin  Az  (2) 

HZ  ~ sin  El  (3) 

These  equations  permit  calculation  of  the  direction  cosines  if  Az  and  El  are  given. 
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TRAJECTORY 

For  u single  shot,  it  is  assumed  that  Az  and  FI  are  uniformly  distributed  between  the 
limits  Azmin,  AzmaX’  F1m1N*  and  FImaX-  FI  is  found  by  selecting  an  ru  (random 
number)  uniformly  distributed  between  zero  and  one  and  using  it  in: 

E1  = E‘{dlN  + ru(E1MAX’BMlN>  (4> 

Az  is  found  by  selecting  anotlter  ru  for  use: 

Az  = AzM1N  + ru(AzMAX‘ ^MIN*  <5> 

The  direction  cosines  jix<  My.  and  can  be  calculated  from  equations  (1),  (2),  and  (3). 

The  Vsh  bas  magnitude  V§H  and  is  expressed  in  equation  (6),  where  i,  j,  and  k 
represent  the  unit  vectors  in  the  X,  Y,  and  Z directions,  respectively.  The  (target 
velocity)  has  magnitude  V-p , is  in  the  X direction,  and  appears  in  equation  (7): 


VSH  = VSHMX  i + VSHPy  j ' vSHMz  k <6> 

VT=VTf  (7) 

— ► 

The  Vr  (relative  velocity)  has  magnitude  Vr,  and  is  defined  and  evaluated  as  follows: 
Vr  = Vsh-Vt  = (VshMx' VT)  ^+VsHMyT’+ VSHMZ  k (8) 

VrWVSH-2VSHvT»‘X+vT  (9) 

Let  A,  B,  and  C be  the  direction  cosines  of  Vr  with  respect  to  the  X,  Y,  and  Z axes. 


. VSH  Mx " VT 

A = — v£ 


B 


C = 


VSH  MY 
Vr 

vSHMz 

Vr 


(10) 

01) 

(12) 


The  Dm  is  the  closest  approach  distance  of  the  trajectory  to  tlie  aim  point  which  in 
this  case,  is  the  center  of  gravity  of  the  target.  The  method  used  to  choose  a value  for  D\f 
depends  on  the  weapon  type. 
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For  the  weapons  with  fragmenting  warheads,  an  a a (aiming  sigma)  is  an  input.  The 
value  of  Dm  is  chosen  randomly  from  a normal  distribution  having  mean  zero  and  standard 
deviation  a\. 

The  solid  shot  weapons  are  characterized  by  a table  of  mil  errois  as  a function  of  target 
velocity.  The  Em  (probable  aiming  error  expressed  in  mils)  can  be  found  for  any  Vx  by 
interpolation  in  this  table.  Let  the  altitude  of  the  target  be  H,  and  the  R§  (slant  range)  is 
given  by : 

RS  = H/sin  (F.1)  (13) 

The  Ep  (probable  aiming  error  expressed  in  units  of  linear  measure)  can  be  calculated  from: 

Ep  = 0.001  Rs  Em  (l4*) 

For  normal  distribution,  the  Ep  is  equal  to  0.675  standard  deviation  units.  Thus,  for  the 
aiming  error. 


OA  = Ep/0.675 


(15) 


Given  the  aiming  error,  the  Dm  defines  the  radius  of  a circle  in  the  plane  perpendicular  to 
the  Vr  and  centered  at  the  aiming  point.  The  following  procedure  chooses  a point  (X0,  Y0, 
ZG)  at  random  on  the  circumference  of  the  circle.  This  point  and  the  Vr  defines  thS 
trajectory  for  a single  shot. 

Select  an  ru  uniformly  distributed  between  zero  and  one.  Let  # = 2jt  ru.  A,  B,  and  C 
are  the  direction  cosines  of  Vr  as  calculated  from  equations  (10),  (11),  and  (12).  If  C ^ 1, 
let: 


COS  # =*“ 


Then, 

IfC  = 0, 
Otherwise, 


sin  # = 


Xq  = |Dm  | 1C  cos#  cos  # + sin  # sin  #) 
Y0  = I Dm  | l-C  cos  # sin  # + sin  # cos  #] 

Zo=  Pm|cos# 

Zq  = (-AX0  - BY0)/C 


(16) 

(17) 

(18) 
(19) 
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If  C = 1,  the  following  equations  apply: 

Xo=PM|cos0  (20) 

Yc  = 1 Dm  J sin  <fi  (21) 

Zo  = 0 (22) 


GEOMETRY  OF  SOLID  SHOT  ENCOUNTER 

For  solid  shot  weapons,  a hit  is  said  to  occur  if  the  trajectory  intersects  any  fuel  tank. 
Consider  the  case  of  one  tank  having  dimensions  LT,  WT,  and  HT,  and  centroid  located  at 
(XcG.  YcGi  and  Zcg)-  For  a given  trajectory  there  are,  at  most,  three  faces  of  the  tank 
through  which  it  is  po&ible  for  the  shell  to  enter.  These  can  be  determined  from  the 
direction  cosines  of  the  Vr  as  shown: 


<0 

=0 

>0 

Front 

No  intercept  with  froru  or  rear 

Rear 

Left 

No  intercept  with  side 

Right 

Top 

No  intercept  with  top  or  bottom 

Bottom 

The  procedure  for  determining  whether  a hit  occurs  is  to  find  the  coordinates  of  the 
points  which  represent  the  intersection  of  the  trajectory  with  those  planes  (taken  from 
Table  2)  in  which  the  faces  of  the  tank  lie.  These  points  are  examined  to  determine  whether 
they  fall  within  the  bounds  which  form  the  faces  of  the  tank.  For  example,  the  planes  which 
contain  the  top  and  bottom  faces  of  the  tank  are  parallel  to  the  XY  plane,  and  their 
equations  are,  respectively: 


Z = ZCG  + HT/2 

(23) 

Z = Zcg*HT/2 

(24) 

The  equation  of  the  trajectory  is: 

X - Xp  _ Y - Yp  _ 7.  - Zp 

a b r 

Let  (Xin,  Yin,  Zin)  be  the  intersection  point  of  the  trajectory  with  whichever  plane,  (23) 
or  (24),  is  encountered  by  the  shell  first.  Then,  using  Table  1 and  equations  (23)  and  (24): 


-JN  = ZCG  + HT/2,  C<0 
Zjn  = ZCG  * HT/2,  CX) 
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Substituting  this  value  of  Zjn  into  equation  (25): 

Xjn  = A/C  (Zin  - ZQ)  + X0 
Y1N  = B/C  (ZiN  - Z0)  + Y0 

If  C is  equal  to  zero,  the  trajectory  is  parallel  to  the  top  and  bottom  faces  and  the  point 
(Xin.  Yin,  Zin)  does  not  exist.  This  point  indicates  a hit  on  the  top  or  bottom  face  of  the 
tank  if  the  following  conditions  are  satisfied : 

Xcg-LT/2<Xin<XCG  + LT/2 
Yeti  - WT/2  < Yin  < Y(j(j  + WT/2 

If  a hit  is  not  found  on  the  top  or  bottom  face,  the  sides  and  the  front/rear  faces  are 
checked  by  a procedure  similar  to  the  above.  If  a valid  entry  point  is  found  for  one  of  the 
three  possible  faces,  the  other  three  faces  are  checked  to  determine  the  exit  point.  Based  on 
the  percent  of  fuel  remaining,  it  is  determined  whether  the  entry  and  exit  points  are  above 
or  below  the  fuel. 


GEOMETRY  OF  F RAGMENTING  WARHEAD  ENCOUNTER 

The  method  of  treating  the  fragmenting  warhead  is  to  determine  the  burst  point  and 
solve,  using  an  iterative  procedure,  the  equation  which  relates  distance  and  time  traveled  for 
the  fragments.  Having  solved  this  equation,  and  knowing  the  fragment  density,  it  is  possible 
to  calculate  the  expected  number  of  hits  on  a tank  and  the  probability  of  a kill. 

The  standard  deviation  of  the  fuzing  error  along  the  trajectory  is  op.  The  aiming  point 
for  this  type  of  weapon  is  not  assumed  to  be  the  center  of  gravity,  but  is  input  as  (XpusE* 
YpuSE*  ani*  Zfuse)-  The  error  along  the  trajectory  due  to  fuzing  is  taken  from  the  point 
of  closest  approach  to  (XpusE'  ^FUSE’  an^  ^FUSE)-  This  error,  Dp,  is  chosen  at  random 
from  a normal  distribution  having  mean  of  zero  and  standard  deviation  ap.  For  a burst 
point  located  by  X*,  Y*,  Z*: 

X*  = X<>  + YpusE  + OF  0*X  ‘ VT/VSH) 

Y*  = Y0  + YpusE  + DFMY 

Z*  = ZG  + ZFUSE  + Dp  mz 

The  explosion  of  a static  fragment  warhead  yields  a characteristic  spectrum  of 
fragment  mass,  angular  density,  and  emission  velocity.  The  explosion  of  a moving  fragment 
warhead  alters  this  spectrum  by  virtue  of  the  velocity  of  the  velocity  of  the  projectile.  It  is 
necessary  to  determine  the  interaction  of  this  altered  spectrum  with  the  target.  The 
relationship  between  speed  and  direction  of  the  projectile,  and  the  speed  and  direction  of  an 
emitted  fragment  are  derived  using  Figure  4.  The  Vp  (fragment  emission  velocity)  and  the 
angles  are  those  observed  in  a static  explosion:  while,  through  vector  addition,  V0 
(observed  fragment  velocity)  and  angle  7 occur  in  a dynamic  explosion. 
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The  law  of  cosines  applied  to  Figure  4 yields  a quadratic  equation  for  VQ: 
vo  " VS11  cos  7 ± y(VsH  cos7)2  + V|-V|H 

IfVSH  < Vp,  as  is  usually  true,  the  negative  root  leads  to  a negative  velocity  which  is  ruled 
out.  Thus,  equation  (26)  is  valid  and  emission  velocity  is  single-valued  for  a given  7 when 
VSH<  Ve: 

Vo=  VSHcos'lr+  \/(VsH  cos  7)2  + VE'VSH  (26) 

Another  relationship  results  from  Figure  4 by  summing  vector  components  in  the  Vsh 
direction: 


cos  0 


V0  cos  7 - Vsh 

Vp; 


(27) 


The  fragment  ballistics  must  be  considered.  A stationary  x,  y,  z coordinate  system  is 
employed.  This  coordinate  system  is  defined  to  coincide  with  the  X,  Y,  Z system  at  the  t=0 
(time  of  the  explosion).  The.  explosion  point  in  the  x,  y,  z system  is  at  (x*,  y*,  z*),  where: 


x*  = X* 
y*  = Y* 
z*  = Z* 
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The  target  is  represented  by  a set  of  fuel  tanks,  and  the  location  of  each  tank  is 
designated  by  the  coordinates  of  its  center  of  gravity.  It  must  be  determined  whether  each 
tank  is  hit  by  the  fragment  spray,  and  if  so,  what  the  expected  number  of  hits  will  be. 

Using  this  procedure  for  one  tank,  consider  a tank  located  at  point  (XcG*  YCG>  ZCG)- 
At  t=0,  the  tank  coordinates  in  the  stationary  system  are  (x,  y,  z)  = (XqG>  YCG>  Zcg)-  The 
tank  is  moving  in  the  +x  direction;  therefore,  the  hit  point  occurs  one  time-of-flight  later  at 
point  (x,  y,  z)  = X|j,  Ycg>  Zcg)-  At  this  point  the  fragment  has  traveled  a distance  L, 
where: 


*H  = XcG  + VT1 

(28) 

L = \/ (XII  • x*)2  + (yCG  * y*)2  + (ZCG  * z*)2 

(29) 

The  direction  cosines,  with  respect  to  the  x,  y,  z axes,  of  the  line  from  explosion  point  to 
hit  point  are  px,  P y,  and  pz: 

0x=  HL 

(30) 

%=  °Ly 

(31) 

(32) 

With  the  direction  cosines  of  each  vector  known,  the  angle  7 (Figure  4) 
taking  the  scalar  product  of  the  vectors: 

given  ,fcy 

cos  7 = Px  Px  + My  Py  + Pz  Pz 

(3>  ■ 

The  distance-time  relationship,  which  describes  the  fragment  travel,  can  be  derived.  By 
equating  the  inertial  and  drag  forces  for  geometrically  similar  bodies,  it  can  readily  be  shown 
that  the  logarithmic  derivative  of  velocity  with  travel  is  proportional  to  the  air  density,  and 
inversely  proportional  to  any  characteristic  length  of  the  body.  The  proportionality 
constant  is  determined  by  the  drag  coefficient,  the  mass  density  of  the  body,  and  the 
geometrical  shape.  Thus,  for  fragments  having  some  characteristic  mass  spectrum,  a sea  level 
slowdown  constant  (k)  may  be  introduced  whose  value  will  be  independent  of  mass: 

dlnV=  k p 

dL  m 1/3 

where  V is  the  velocity  of  a fragment  after  traveling  a distance  L,  p is  the  relative  air  density, 
and  m is  the  fragment  mass.  This  may  be  integrated  at  constant  drag  coefficient  to  yield  the 
velocity-distance  equation  for  fragments,  where  V0  is  the  initial  fragment  speed: 


In  (V/V0)  = - k p L/m1/3 
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This  equation  is  integrated  once  more  to  obtain  the  desired  result: 
B V0t  = eBI~l 


(34) 


where 

B = k p/m  1 /3 

The  criterion  for  a hit  and  all  corresponding  properties  are  determined  by  simultaneous 
solution  of  equations  (26),  (28),  (29),  (30),  (31),  (32),  and  (34).  No  analytic  solution  to 
this  system  has  been  found,  but  an  iterative  numerical  solution  can  be  employed. 

The  numerical  method  is  the  Newton-Raphson  technique.  The  formulation  of  this 
method,  as  applied  to  fragments,  has  been  tested  and  found  to  give  rapid  convergence  even 
with  inputs  which  were  known  to  be  troublesome  by  previous  methods.  The  Newton- 
Raphson  method  obtains  the  roots  of  F(t)  = 0. 

The  classical  method  takes  the  Jth  estimate  of  the  root  t-J  and  extracts  the  (J+l) 
estimate  t(^+^  by  means  of  equation  (35): 

t(J+l)  = t(J)-  (35) 

F(tJ) 

The  procedure  is  repeated  until  successive  estimates  are  considered  to  differ  negligibly. 
Application  of  this  method  to  the  fragment  ballistics  required  suitable  choice  of  F(t).  The 
travel-time  relation  (34)  is  chosen  and  rewritten  as  the  F function: 


F(t)  - In  ( 1 + B VQ  t)  - B L 
Differentiation  yields  the  time  derivative: 


F(t)  = B 


Vp  + Vot  ■ 
1 + B V0  t 


(36) 


(37) 


The  two  time  derivatives  in  equation  (37)  must  be  evaluated  from  the  other  relations  that 
must  be  satisfied. 


Differentiation  of  equation  (26)  at  constant  Vg  yields: 


V0  = 


-Vo  VSH  sin  7 
V0  - VSH  cos  7 


7 


(38) 


Differentiation  of  equation  (28)  yields: 
xh  = Vj 


(39) 


Differentiation  of  equation  (29),  and  use  of  (30)  and  (39)  yields: 

L = 0X  vT 


(40) 
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It  can  be  seen  that  while  L is  known  in  terms  of  non-derivatives  V0  is  known  in  terms  of  7. 
Thus,  to  complete  the  evaluation  of  equation  (37),  it  is  necessary  to  calculate  7. 


Differentiation  of  equations  (30),  (3 1 ),  and  (32),  and  using  (39)  yields: 

.*  VT-0xL 
Px  “ T 


• -jjy  L 


(41) 

(42) 


(43) 


Differentiation  of  equation  (33)  and  employing  (41),  (42),  and  (43)  gives: 

. L cos  7 - px  Vj 
**  L sin  7 


(44) 


• • 

The  final  results  are  obtained  by  eliminating  L and  7 between  equations  (37),  (38),  (40), 
and  (44): 


F(t)  = B 


vo  VSH  VT 

Mx-Pxcos  7 ' 

L 

vo‘  VSH  cos  7 

(45) 

(46) 


Thus,  the  Newton-Raphson  method  for  fragment  ballistics  employs  the  system  of 
equations  (35),  (36),  (45),  and  (46). 

The  method  used  to  determine  the  hit  point  of  a given  fuel  tank  for  a fragment  of  a 
given  mass  is: 

a.  Estimate  time-of-flight  (t). 

b.  Calculate  hit  point  coordinate  (xh)  from  equation  (2S). 

c.  Calculate  fragment  travel  (L)  from  equation  (29). 

d.  Calculate  direction  cosines  (0X,  /3y,  and  (3Z)  of  fragment  velocity  from  equations 
(30),  (31),  and  (32). 

e.  Calculate  dynamic  fragment  emission  angle  (7)  from  equation  (33). 

f.  Calculate  dynamic  fragment  emission  velocity  (V0)  from  equation  (26). 
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g.  Calculate  F(t)  from  (26)  and  V0  from  equation  (46). 

h.  Calculate  F(t)  from  equation  (45).  If  F(t)  < 0,  there  is  no  hit  point  and  the 
fragment  has  missed  the  fuel  tank. 

i.  Make  new  estimate  of  time-of-flight  (t)  using  equation  (35).  If  the  new  t<0,  score  a 
miss;  otherwise,  compare  with  the  previous  value  of  t.  If  two  successive  values  are  in 
agreement  (e.g,,  result  in  a difference  of  less  than  0.5  foot  in  the  hit  point),  the  process  is 
considered  to  have  converged  and  the  most  recent  value  of  t is  saved  as  the  solution. 
Otherwise,  iterate  again  at  step  c. 

It  is  important  to  make  a good  estimate  of  the  time-of-flight  for  step  a.  The  method  used  for 
the  first  estimate  is  to  take  the  analytic  solution  for  the  case  of  zero  fragment  slowdown. 

The  distance-time  relationsliip  for  the  zero  slowdown  case  is; 

L = VQt  (47) 

With  a substantial  amount  of  algebraic  manipulation,  equation  (48)  is  combined  with  (26), 
(28),  (29),  (30),  (31),  (32).  and  (33)  to  yield. 

Kjt2  + K2t  + K3  = 0 (48) 

where 

Kl  - Vr  - V£ 

K2=-2  Vsh  [(#ix-  *xCC*x*)+My  <YCG  * y*)  + H (ZCG  - Z*)J 
K3  = (XCG  - x*)2  + ( Ycg  - y*)2  + <ZCG  * z*)2 


The  solution  to  this  equation  is: 


t = 


*2_ 

2Ki 


The  smaller  positive  value  of  t from  the  above  solution  is  used  as  the  initial  estimate  for  the 
iterative  procedure. 


When  the  iterations  have  been  found  to  converge  for  a particular  fuel  tank,  the  value  of 
cos  0 is  determined  from  equation  (27).  This  value  is  compared  with  the  limits  of  the  static- 
fragment  spray.  If  the  value  lies  outside  these  limits,  the  shot  is  scored  as  a miss.  If  cos  0 lies 
within  the  bounds  of  the  static  fragment  spray  limits,  it  is  necessary  to  calculate  the 
fragment  density  resulting  from  the  dynamic  explosion. 
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The  model  assumes  that  the  fragment  density  resulting  from  a static  explosion  is 
uniform  between  the  fragment  spray  limits.  Let  £ be  this  density  expressed  in  fragments  per 
steradion.  As  a result  of  the  rotation  of  velocity  vectors  due  to  the  projectile  motion,  the 
fragment  density  in  the  static  case  at  angle  0 is  not  equal  to  the  dynamic  density  at  angle  7. 
This  effect  may  be  calculated  using  the  geometry  of  Figure  4. 


The  solid  angle  subtended  by  the  conical  shell  between  0 and  0 + d0  is  dc oq\ 

(Vfcdtfl  [2*  Vg  sin  0*1  . 

d U3Q  =* Jss  2t  sin  0 d 0 

V£ 

Similarly,  the  solid  angle  of  the  conical  shell  between  7 and  7+d  is  dw<y‘. 


dui-i 


[Vo  <b]  O V0  sin  7]  , . 

= It  sin  7 d7 


V2 

vo 


Assuming  that  adjacent  rays  satisfy  the  geometry  of  Figure  1,  the  number  of  fragments 
in  each  of  these  conical  shells  is  the  same;  thus,  E is  a measure  of  the  change  in  fragment 
density,  where: 


£ _ s*n  7 d 7 

du I#  f 1 0 d 0 


(49) 


*DYN  ■ i'E 

To  calculate  the  value  of  E must  be  derived  for  use  in  equation  (50). 

The  speed  ratio  G is  defined  as: 


(50) 


G = 


VE 

VSH 


Elimination  of  VQ  between  equations  (26)  and  (27)  gives: 


1 . 1 , _ G cos  0 + 1 

cosz  7 + G - * 1 * — — - — • cos  7 

cos  7 

Squaring  and  rearranging  terms  yields: 

( 1 + G cos  0 )2 


cos2  7 = 


(G^  + 2 G cos  0 + I ) 
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Differentiating  equation  (5 1 ) results  in;  ; ; ^ 


dy  _ G-  sin  6 (\  + G cos  d)  (ti  + cos  (?) 

.. . dO  ■'  cos  7 sin  7 (G2  + 2 G cos  0 + ' w 


Substituting  this  result  in  equation  (40): 


• j *»'  ' 


G^fl  + G cos 0) (G + cos  0) 
E •=  * — 


V 

. V* 


jv.  . C06.7  (G2  + 2G  cos  0 4 i')2 


Ojas  y between  equations  (51)  and  (52).  Based 


E = 


G2|C  + cos  6\ 


(G2  + 2 G cos  $ + 1)3/2 


The  fuel  tank  is  treated  as  a rectangular  parallelepiped:  therefore,  there  are  a maximum 
of  three  faces  that  can  be  hit  due  to  one  explosion,  identification  of  the  three  faces  can  be 
achieved  by  the  use  of  relative  velocities. 


VH1T  (fragment  speed  at  the  time  of  the  hit)  is  found  to  be: 
^HIT  = ^oe 


The  striking  fragment,  target,  and  relative  Vjq£T  velocities  can  now  be  expressed: 


VHiT  " VHIT  2x  ‘ + VHJT  % J + VHIT  h k 

* vp 

VNET  = ^HIT  ' Vt  = ( VHIT  Px  * Vj)T'+  Vhit  Py)  + VHit  Pi* 
The  VNET  (net  striking  speed)  is  the  magnitude  of  the  velocity  VjqET: 


VNET  =yjv&l T - 2 Vj  Vhit  Px  + v| 
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1 %e  signs  of  the  Vj^j:  j components  identify  the  struck  aspects  as  follows: 


% 


<0 

=0 

>0 

0Z 

Top 

No  strikes  on  top  or  bottom 

Bottom 

l£ft 

No  strikes  on  sides 

Right 

Yt 

Front 

No  strikes  on  front  or  rear 

Rear 

Vhit 


The  number  of  fragment  hits  on  each  of  the  three  aspects  remains  to  be  calculated. 
Calculation  of  the  fragment  density  in  target  coordinates  appears  tedious  and  possibly 
difficult.  Therefore,  the  approximation  is  made  that  the  number  of  hits  can  be  calculated  on 
a static  target.  This  should  be  an  excellent  approximation  if  VHIT  » Vj!  and  probably  not 
too  bad  for  most  cases  to  be  encountered.  Using  the  static  target  concept,  Table  2 is 
modified,  replacing  Vt/Vhjj  by  zero. 

Consider  a static  target  placed  in  a constant  density  pulsed  beam  of  particles  emitted 
from  a static  point  source.  To  be  definite,  consider  the  right-left  aspect  only.  The  actual  area 
of  tike  aspect  is  Ay,  while  the  area  component  normal  to  the  beam  is  Ay|0y|.  The  solid  angle 
viewed  by  the  point  source  is  approximately  Ay|/?y  |/L2.  If  the  separation  L is  quite  small, 
this  will  give  a large  overestimate  of  the  solid  angle,  but  this  is  not  important  since  kill  will 
be  achieved  for  small  L.  Thus,  the  number  of  hits  on  the  aspect  is  £dYN  Ay|j3y|/lA  From 
this  and  equation  (50),  the  following  result  is  generalized  for  the  i^  aspect: 


SAMI 

"i  = EL2 

where  nj  is  the  number  of  hits  on  the  aspect. 


(54) 


If  N is  the  total  number  of  fragments  emitted  by  the  warhead,  the  probability  that  a 
fragment  selected  at  random  scores  a hit  on  the  i^1  aspect  is: 

Pi=^- 
1 N 

Considering  each  fragment  to  be  independent  of  other  fragments  in  hitting  a fuel  tank,  the 
PH  based  on  the  appropriate  three  faces  can  be  formulated  by: 

1-PH  = n(l-pi)N  = n(l-Pi)('  pP  ('ni) 

three 

faces 


[\ 

>■ 

~I1  e *n>  = 

t 

* 

three 

k- 

* • 

faces 

i 

t 

i 

! 

: 
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where  the  approximation  is  good  only  for  pj  « 1.  Substituting  from  equation  (54),  the 
final  form  is: 


Pjl  = 1 - e ^ EL2 


three 

faces 


(55) 


Similarly,  the  Pgp  (probability  of  a hit  below  the  fuel  level)  and  the  PAp  (probability  of  a 
hit  above  the  fuel  level)  can  be  calculated  from: 


PBF  = I - e 


v iAjifti 
^Ei T 


PAp=  1 -e 


-2 


gAjJftl 

EL2 


(56) 


(57) 

If 


where  a(  represents  that  portion  of  the  itfl  aspect  area  which  is  below  the  fuel  level,  and  Aj* 
represents  that  portion  of  the  i(h  aspect  area  which  is  above  the  fuel  level. 

The  probability  of  a leak  is  considered  to  be  equal  to  the  Pgp.  The  PEp  (probability 
that  an  external  fire  occurs  given  that  a leak  exists)  is  calculated  from: 

PEF  = PbfDEF  (58) 

Dpp  is  a degradation  factor  which  is  dependent  upon  the  altitude  at  which  the  encounter 
takes  place  and  the  temperature  of  the  fuel  in  the  tank: 


DEF  = 


0.3 

0.3(T/45) 

0.3(1. 2-.00002H) 
>0.30/45)  (1.2-.00002H) 


H>60,000,  or  H<1 0,000  and  T<0,  or  T>45 
H>60,000,  or  H<1 0,000  and  0<T<45 
1 0,0Q0<H<60,000  and  T<0,  or  T>45 
1 0,0G0<H<60,000  and  0<T<45 


where  T is  the  fuel  temperature  in  degrees  Fahrenheit,  and  H is  the  altitude  in  feet.  This 
relationship  for  the  degradation  factor  is  based  on  limited  data  for  wet  hit  test  results2 

The  Ppp  (probability  of  an  internal  fire/explosion)  is  considered  to  be  zero  unless  the 
fuel-to-air  ratio  in  the  ullage  space  is  within  the  flammability  limits  for  the  particular  fuel 
being  used.  If  the  fuel-to-air  ratio  lies  within  the  flammability  limits  (e.g.,  0.013  to  0.08  for 
JP-4),  this  probability  is  given  by: 


pFE=  PAFDFE 


(59) 


where 


H<  10,000 


^ (.00000769m1/2  VNET 

Dfe  " j / \ 

( .00000769m1/2  VNET  \2.5e -■00092H)  H>10,000 
This  relationship  was  derived  by  fitting  curves  to  data  supplied  by  BRL. 

Ballistic  Research  Laboratory.  Fragment  Firings  Against  Aircraft  Fuels  at  Simulated  Altitude,  by  W.R,  Harris. 
Aberdeen  Proving  Ground,  MD,  BRL,  October  1953.  (BRL  TN  828,  publication  UNCLASSIFIED.) 
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The  last  type  of  damage  mechanism  which  is  considered  by  the  model  is  damage  due  to 
hydraulic  ram.  This  mechanism  is  treated  from  an  energy  density  standpoint.  If  the  energy 
density  of  the  fragment  spray  on  a particular  fuel  tank  is  above  a threshold  value,  ram  is  said 
to  occur  on  this  tank  for  tliis  explosion.  The  energy  density  of  the  spray  is  calculated  from 
the  relationship: 


eRAM  = 


mV$ET 

€ ' 

2 

.EL2 

APPLICATIONS 

The  model  can  be  used  to  examine  variations  in  the  type  and  intensity  of  threat  as  the 
mission  profile  changes.  The  vulnerability  characteristics  of  the  aircraft  vary  with  time, 
maneuver  history,  threat,  and  threat  exposure.  For  example,  aircraft  which  penetrate  and/or 
deliver  ordnance  at  low  altitude  may  be  exposed  to  a greater  variety  of  hostile  weapons  than 
high  altitude  bombers.  Any  given  weapon  system  may  be  exposed  to  a wide  variation  in 
lethal  threat  as  its  flight  profile  is  changed.  Exercise  of  the  model  will  reveal  the  relative 
severity  of  the  threats  and  indicate  potential  phases  for  laboratory  testing. 

In  the  case  of  a large  bomber,  for  example,  it  is  not  immediately  obvious  whether 
concern  should  be  directed  at  air-to-air  missiles  in  the  cruise-out  phase,  or  at  light  AAA 
weapons  in  the  low  altitude  approach.  At  high  altitude,  ullage  spaces  tend  to  be  oxygen- 
poor  and  quite  cool,  which  inhibits  propagating  fires.  At  low  altitude,  the  fuel-air  mix, 
particularly  immediately  following  descent,  may  reach  near-optimum  flammability,  and 
even  an  otherwise  minor  threat  may  become  lethal.  Aerodynamic  heating  late  in  a low 
altitude  phase  may  produce  flammable  mixtures. 


CONCLUSIONS 


The  vulnerability  model  presents  a system  for  studying  the  dynamic  interaction 
between  fuel  state  and  hostile  threat.  Previous  systems  have  studied  the  fuel  system  statisti- 
cally, with  dynamic  treatment  of  wec^ons  only. 

Theie  are  several  limitations  in  this  model.  Ullage  spaces  were  assumed  to  be  homoge- 
neous. The  effects  of  vent  geometry,  slosh,  vibration,  and  splash  caused  by  impact  were  not 
treated.  There  was  no  integration  of  the  fuel  system  into  the  aircraft  structure  (masking  and 
shielding  by  other  components).  Secondary  ignition  sources  were  only  crudely  treated.  Tank 
geometry  was  limited  to  rectangular  shapes.  Round  breakup  and  ricochet  were  not  treated. 

This  project  was  exploratory  in  nature.  The  results  achieved  represent  an  improvement 
in  treating  fuel  system  vulnerability.  Vulnerability  can  now  be  calculated  as  a function  of 
the  mission  style,  as  opposed  to  the  single  point  computations  previously  possible.  This 
represents  a large  increase  in  the  realism  of  vulnerability  computations. 
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Appendix 

PROGRAM  LISTING  AND  SAMPLE  CASE 
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WSFT  PROGRAM 


UOBUSFT  «CM50000,1160,I01QO«P6. 

FTN(OPTaQ) 

LGO( « < « » TNK 1 ) 

AT ALOG ( TNK 1 , TNK1 , RP=1 U • KN=1 ) 

9 END  OF  KEtOKO 

PROGRAM  MAIN ( INPUT (OUTPUT . TAPES* INPUT . T APE6*0UTPUT .TAPE? > 
COMMON/T  ITU/C  10  (12 ) 

COMMON/S  I IK/  RGAS.tHWA.LMWF .CP A. CPF  «HJ(3) .UU.6HtUEWXD.ULHT. 
lOECHF «ZOG ,AV  «CP ,EMW  «UV .CUELT A. K TANK, GALL  «T VENT 
COMMON /EULIN/  DUMMY  ( 9 ) . EM - T A ,Z ,EHV . tHOG .EHEV . EMCD . V . GAL 
C0MM0N/0UT6 AS/BET  A .CONI . C0N2 .HHOLIQ » EMUISQ  « VJMMOO , EHQOOT 
X.VLIO.tMuIS.EMUISE 

COMMON/HMITS/TO.THAX.uT.THRNT.DTPRNI .TIME.M16 
COMMON/ T ABLS/TaB <100.2.61.  NTAB(6> ,NSAVE(6» 

COMMON  /FILE7/  PAR (250 ) .FLEFT (250) .FTEMP (250 ) .HT (260 J . VT (250) . 
1IIPRNT .GALX.TINC 
DIMENSION  HUH ( 6 ) 

DIMENSION  TF(1U0.2) .TSIDE <100,2 >,TTOP (100.2) ,OUM(l) 

S .ALT  (100.2)  « GALUOT  ( 100 .2  > • PVAPU00.2),  EMINFdOO.2) 

EQUIVALENCE  (TAB(1,1«2)  • TF(l.l)  > 

t .(TAB (1.1. A)  , 1 SIDE (1.1)  ) 

S .(TAB (1.1.4)  . TTOP(l.l)  ) 

$ .(TAB (1.1.5)  . ALT (1.1)  ) 

S .(TAB (1.1. 6)  • GALUOT (1.1)  ) 

$ , (TABU, 1*7)  « PVAP(l.l)  ) 

S .(TABU, 1,6)  . EMINF  (1.1)  ) 

S , ( DUM ( 1 ) « TAB (1.1,1)  ) 

NAMELIST/OAT a/RGAS.EMuA.EMrFvCPA. CPF. TA.HJ.ULLGH.ULWID. 

$ ULHT.OELHF, ZOO. AV.OV.COELTA.KTANK.GALO, BETA, CONI, CON2 

S , 1 WEN  t 

S .TF.TSIUE.TTOP, ALT, GALOOT, PVAP.EM1NF 
S , TO , TMAX ,UT .DTPRNT 

DATA  HDH/1H  , 2HTF , 5HTS1UL .4HTT0P,  3HALT  « 6H6ALD0T « 4HPVAP , 

1 5HEMINF  / 

00  5 1=1,1600 
DUM(1I=  -992.1 
5 CONTINUE 
IIPKNTsQ 
IPKT=0 

50  HEAD (5,900)010 
900  (OHMAT ( 12A6 ) 

IF (EOF (5) 11000.10 
10  KEAU ( 5. DAT  A ) 

TINCsUTPKNI 

GALOsGALO 

TXME=TU 


TPKN  f = TO  +UTP8NT 
DT=  AM1N1 ( OT .DIPRN I ) 

TOL  = AMAX1C  .0001*UT,1.0E-6) 

WHITE (6,910 ) 

910  (OHMAT(1H1,20X* 12H INPUT  TABLES  ) 

UO  7b  1=2.6 
WKITL{6,920)  HDR(I) 

920  (OKMAI (1H0.25X.10HT  ABLE  (OK  . A6/) 

NSAVE ( I )=1 
DO  74  0=1,100 

IF ( TAB (0,1,1)  .EQ. -992.1  .ANO.  T AB (0, 2 . I ) .EO ,-992. 1 (GO  TO  75 
WRITE (6. 930)  J , TAB ( 0, 1 , 1 ) . TAB (0 ,2 , 1 ) 

930  (OKMAT(15X,I5»2L17.7> 

NTAB( I )=0 


O 
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SUBROUTINE  LINK X . Y . N < AHG  < YANG  < NSA VL I 
C 

C LINEAR  INTERPOLATION  ROUTINE 

C 

dimension  mih  y < i » 

NN=N 

XVsAKG 

IF  (XV.GE.X(NN) 160  TO  40 

ifuv.le.xuiiso  to  so 

c 

J=  NSAVE 

IF  I J .LT.  1 .OK.  U .ttT . NN)  J=1 
K=  SlGNU.U«(XV-*|tJm 
S 0=J+K 

IF  I tXV-X(J)  I*  FLOAT (K)  1 10.30.5 

io  iF(K.e.a.-iusj4i 

• Isd-i 

c 

C INTERPOLATION  CALC 

C 

H=X(J)-X< II 
ux=xv-x< i ) 
oy=y<j>-y<i> 

YAR6=  Y ( I ) ♦ UX*QY/H 
NSAVE.  = 1 

return 
30  YAPGsYIJ) 

NSAVEsd 

Return 

40  VAR6=Y(NN) 

RETURN 
50  YAHGsY ( 1 ) 

RETURN 

ENO 
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LINS  9 
LINX  10 
LINX  11 
LINX  12 
LINX  13 
LINX  14 
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LXNI  20 
LINI  21 
LINI  22 
LINI  23 
LINI  24 
LINX  25 
LINI  2t 
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IF (ABC (DEI ) .EU.  Q.I  60  TO  SO 

LESK0062 

C 

COMPUTE  JE1ERHINANT 

LESK0063 

LiU  40  1 = 1, N 

LESK0064 

‘40 

UET  = OCT«A ( I « I ) *S( 1 ) 

LESK0065 

c 

HACK  SUBSTITUTE 

LESK0066 

so 

A1  = A ( N , N ) 

LESK0067 

00  70  J=lgM 

LESK006A 

JJ  s M+J 

LESK0069 

• 

X ( N , J ) = AlN,JJ)/Al 

LESK0070 

IF(N  .EU.  1)  GO  TO  70 

LESK0071 

UO  65  1=2 ,N 

LESK0072 

K = MN-1 

LESK0073 

A2  = A t K t K > 

LESK0074 

IF ( AbC ( A2 ) .LE.  1.0-10)  GO  TO  110 

LESK0075 

B = A (K , JJ ) 

LESK0076 

LL  = K*1 

LESK0077 

00  60  L=LL,N 

LESK0076 

60 

B = 6-A(K,L)*X(L«J) 

LESK0079 

65 

X(K.Jl  = B/A2 

LESK0030 

70 

CONTINUE 

LESK0061 

LA  = 0 

LESK0082 

return 

LESK0083 

100 

LA  = 1 

LESK0084 

Return 

LESK0085 

110 

LA  = -1 

LESK0066 

RETURN 

LESK0087 

END 

LLSK0088 

REAL  FUNCTION  LIN(AftG,N) 

COFiMUN/TABLS/TAtH  100 ,2,6)  , NT  AB  < 8 1 « NSAVE  ( 6 ) 

CALL  LI.'UUAb(l.l.N)  .TABU, 2, N)  «NTAB<N) , AR6. T .NSAVE (N > ) 

LIN=  r 

RETURN 

END 
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c 


b 

10 

c 


15 

C 


16 


20 

22 


20 

29 

30 


SUBROUTINE  LESK (A'X,S.N1,M1.N1X<DET «LA ) 

LESK0002 

DIMENSION  A(NlX.l) .X(NlX.l)  ,S(N1X> 

LESK0003 

THIS  ROUTINE  lb  A SINGLE  PRECISION  LINEAR 

EQUATION 

LESK0004 

SOLVER.  IN  OKUER  TO  CONVERT  TO  COMPLEX  OR 

DOUBLE- 

LCSK0Q0S 

PRECISION,  REMOVE  THE  FOLLOWING 

CARD. 

LESK0006 

ABC  (OEM  = ASS  (OCT) 

LESK0007 

THEN,  IF  DOUBLE-PRECISION  IS  DESIRED,  REMOVE 

LESK0009 

THE  COL  1 C IN  EACH  OF  1 HE  NEXT 

TWO  CAROS 

• 

LESK0009 

DOUBLE  PRECISION  A , X ,S *UET . AS . A1 . A2 ,K . B 

LESK0010 

AUC(UET)  = DABS ( UET ) 

LESK0011 

IF  INSTEAD  COMPLEX  IS  DESIRED,  REMOVE  THE 

COL  1 C 

LESK0012 

IN  EACH  OF  THE  NEXT  TWO  CAROS. 

LESK0013 

COMPLEX  A , X ,0£T ,A1,A2,K,B 

LCSK0014 

ABC(DET ) a CABS ( DET ) 

LESK0015 

u = Wl 

LESK0016 

h a Ml 

LESK0017 

MN  = N+l 

LESKOOIB 

NM  = N-frt 

LESK00I9 

GET  SCALE  FACTORS 

LESK0020 

DO  10  1=1, W 

LESK0021 

S(l>  = ABC (A(I,1)) 

LESK0022 

DO  5 J=1 , N 

LESK0023 

DA  a ABC ( A ( I , J ) ) 

LESK0024 

IF ( UA  ,LE,  S(I)>  GO  TO  5 

LESK0025 

S(I)  a DA 

LESK0026 

CONTINUE 

LESK0027 

IF(SU)  ,EU«  0.)  GO  TO  100 

LESK002B 

CONTINUE 

LESK0029 

SCALE  ROWS 

LESK0030 

DU  15  1=1, N 

LESK0031 

AS  = l./MD 

LESK0032 

DO  15  d=l,NM 

LESK0033 

A ( I , J ) = AS* A ( I , J ) 

LESK0034 

START  TRIANGULAR  1,2 A T ION  PROCESS 

LESK0035 

IF (N  ,EQ.  1)  GO  TO  35 

LESK003G 

NO  = N-  1 

LESK0037 

DU  30  1=1, NO 

LESK0038 

K = 1 

LESK0039 

DA  a ABC  (Ad,  II) 

LESK0040 

UO  16  J=1,N 

LESK0041 

OB  a ABC  ( A ( <J  « I ) > 

LESK0042 

IF ( UD  *LE,  OA)  GO  TO  16 

LESK0043 

Kao 

LESK0044 

UA  a OB 

LESK0045 

CONTINUE 

LESK0046 

IF ( UA  ,EO.  0.)  GO  TO  30 

LESK0047 

IF ( K .EU.  1)  GU  TO  22 

LESK0046 

DO  20  Jal,NM 

LESK0049 

B a A ( K , J ) 

LESK0050 

A(K,J)  a A ( 1 , U ) 

LESK0051 

A ( I , J ) a B 

LESK0052 

DET  a -DET 

LESK0053 

11  = 1*1 

LESK0054 

DO  29  U=II,N 

LESK005S 

R = A(J. 1) /A (I, I) 

LESK0056 

□ U 26  K = I I * NH 

LESK0057 

A ( J , K ) a AMiK)-K*A(I,K) 

LESK005B 

CONTINUE 

LESK0059 

CONTINUE 

LESK0060 

IF  t ABC  t A ( N *N ) ) ,LE.  1,0-10)  GO 

TO  110 

LESK0061 

i • iL 


•MaaiiaWMlii 
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SUitf.OUl  lrjt  K41T1L 

CCf  HOU/LlHlTS/IO.TrtAX.UT  .TPHWT iOTPKNI , UME.M16 
COC'HGU/STIX/  KGaS.LMWA.EMwF  .CPa.CPF  .hJt  3 1 .ULL6H.ULWI0  *UCHT  * 
IDELHF.ZOG.AV.CP.EhW.UV  , tutu  T AilU  ANA  « GAL  0 » TVEWT 
C04PUINI/EUL1N/  QUHHT  t9>  .EK . T A . Z « tMV  .EH0l»  . ERE V .ERCO  . V .GAL 
CORRON/OUt&AS/bETA.CONl.LOhZ.RHOHa.iRUISO.SURHDO.EROt'OT 
I.VLIU.ENQIS.LRUISE 
PL  At.  LlN 

NAMELIST/CMNPT/  RGAS  <ERmA .EHtaF .CPA. CPF .TA.HO. 

1 UULOH.ULtalU.ULHT.OLLHF « ZOG . AW . UW . CUELT A .KT ANK . GALG 


SUPMUO  = U.O 
CftV  =0. 

EflLV  = 0. 

EP-iOG  = 0. 

tncti  = u. 

GAL  s 0. 

TA  = TA  *459.7 
IF  ( AT AimK  .EU.  1)  GO  TO  100 
V = DLL &M * JL» I C * ULKT 
GO  Tt!  2UU 

100  V = 3,14159265*  ULWlO**2*ULLGH/4.  -GAL0*Z31 . / 1728 . 

200  ALT=LiN(T 1ME.51 
ALT  = ttLT*1000.0 

CALL  ATAiOSt  ALT  .TALT  ,PR.0UMPnUuni.UtJf12> 

TFS  LlNtTIht.2) 

PVAP  = L1N<  TF.7) 

TF  s TF  *459.7 
PPFULL  - H4AP*144, 

PPA1N  s PH  - PPFULL 

RHU  s (PPFUEL*EMWP  * PPAlK*Eh4A)/lRGAS*TA> 

Eh  s HHU*V 

2 = PPAIK*EMWA/(RHO*RGAS*TA) 

VLIQU  = GAL 0*231. /1728. 

EhWOG  = ERWA 

ErtUISu  = (BE  T A *EMWOG*VLi WO *PP AIK )/( ,797*453, *2116.224 1 


C 

C 


WRITE  (fc.CKifJPT) 

RETURN 

END 
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SUbKuUTIWL  tULbH 

COMMUN/Ll'UrS/IO.TrtAX.UT.TPRNl .IMPRNI , tXME.Hlfc 
OIMCNSIOM  X(9ltXU0Uy> 

COMION/tULin/  XUOT < X 

CGiSilO.J/OUTGAS/ULT  A.CUNl.LOU2.RHOLiQ«EMUISO.SUMHOO*EMOOOT 
l.VLXQ.£MUIi>.EMUli>E 


OU  1U  1 = i.9 

x<n=  x ( i.  > ♦ xuut  < x ) «dt 
continue 

SUnrtbO  s SUMMOO  ♦ EMDUOT*  Ul 

RETURN 

ENl) 


SUbKuUT INE  FUNCT 

CUMNUN/ST  IR/  HGAS,EMGA,£MwF*CPa.CPF«HJ<3) ,ULLGH,ULGID.ULHT. 
lOLLHF  i20GtAV'CP«EM«iUV  tCUELT  A «M  T ANR  t GALO  « TVENT 


COMMON/EULXN/  MOOT . TOOT .200T . MOOTV.  MOOTOCt  HDOTEV.  HUOTCD.VOOT- 
1GALOOT. 

2EM, I A »2 1 EMV  »EMGG  «EMEV iEMCO • V » GAL 
COMMON/ETC/  ALT , AJ13) ,TJ< 3) .RHOV.PR.PPAIR.PPFUEL.PPFLQS. 

1 AMAtH,UALTUT,PUOT,OF , MASR AT « UVENT 
COMM&N/OUTGAS/bETA»CONltCON2» KHOL1Q « EMUISO • SUMMOO tEHQDOT 
l.VLIG.EMjXS.ERUlSE 

COKMGN/LIfUTS/TO.TriAX,l)T.TPRNT,OTPRNT,  TIME*  Rib 
REAL  MO . T «MOOTV ,MDOTEV . MOOTED .HUOTOG , MASR AT 

NAMELIST/ANS/  TIMEtALT.  MOOT , TOOT .ZUOT . MOOTV.  MOOTOG,  HDOTEV. 

1 MOOT EU.VDOT. GALOOT ,EM , TA , 2 . EMV .EMQG .EME V ,EHCD . V » GAL .PR . POOT . 
2PPAIK,PPFuEL.PPFLQS.UVENT,OF«MA5RAT 


CALL  EULER 
TIME=TXME^OT 
CALL  DERXV 


RETURN 

ENU 
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ELZ  a ELH*G/60 

OHOZ  = 0.0 

EH  a WHO 

C CHECK  THS  SLOPE  AND  CALCULATE  PRESSURE 

IF  t ELH  .EU.  0.0)  00  TO  S 

C NON  - ZERO  SLOPt  PKESSURE  EQUATION 

A ( 4 I = PHt J)*lTH| J)/THSI**(OHHS/ELH) 

GO  TO  S» 

C ZEHU  SLOPE  PRESSURE  EQUATION 

5 A ( 4 I = PH  ( J ) *EXP I 6HRS* (HG  ( J ) -H) /THS ) 

GO  TO  9 

C THS  LINEAR  WITH  2.  SEARCH  MATRIX 

6 00  7 I a 2 ,.14 

J = I ♦ 8 

K a I - 1 

IF  (ZHII)  .6E.  2)  GO  TO  0 

7 CONTINUE 

C CALCULATE  THS,  SLOPE,  ANO  STUFF 

8 ELZ  a (TMIJ+1)  > TH  ( J ) ) / ( ZM I K> 1 ) - ZH|K)) 

THS  a TH(J)  ♦ ELZ*  I 2 - ZHOU) 

DHUZ  a ( WH I K ♦ 1 ) - WH(K) )/(ZHIK*l>  - ZH(K)) 

EH  = Wh(K)  ♦ UM0Z*tZ  - ZH(K)I 
ZLZ  a Z - THS/ELZ 

C PRESSURE  EQUATION  FOR  THS  LINEAR  WITH  2 

A I 4 ) a PH<U)*EXP(GHHS/ELZ*IKU/(K0*2LZ> )**2*( (Z-ZHlK) )• 

1 (K0,ZL2)/|RU*Z)/IH0+ZH(K))  - ALOG ( THS* ( RU+ZH ( K ) 

2 »/TH( J)/-(R0*Z) ) ) ) 

C CALCULATE  SOUNU  SPEEU  ANO  DERIVATIVE 

S All)  a 49,0221G4*SUKT (THS) 

A 1 2 I a 0 »5*ELZ/THS 

C CALCULATE  UENSITY , DERIVATIVE,  ANO  PKESSURE  UEHIVATIVE 
A(b)  a GHRS* A ( 4 ) /GO /THS 
A ( 7 ) a - ( A ( 6 ) *G/A ( 4 ) * ELZ/TMS I 
Alb)  a - A I G I *G 

C CALCULATE  TEHPERATURf . DERIVATIVE,  AND  LEAVE 

A ( 6 ) a EM*THS/*H0 
AID)  a <EH*ELZ  ♦ THS*0HQZ ) /WHO 
10  At)  a A ( 8 ) 

A4  a A 1 4 ) 

A 1 a All) 

AG  a A I G ) 

AS=AIS) 

RETURN 

END 
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r 


t 


» 


•* 


* 


« 


» 


i 


c 

c 

c 

c 

c 

c 

c 

c. 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SOPH OUT  1 NE.  A I MO  S (A4.Afi.A4.Al.Ab.A5) 

THIS  ROUT InL  CALCULATFS  ATMOSPHERIC  PROPER! IES  OF  THE 
US  STANUARU  ATMOSPHERE , 1962 , ASSUMING  AN  INVERSE  SOUARE 
GRAVITATIONAL  FItLU.  THIS  ASSUMPTION  YIELDS  DATA  THAT 
AGREES  WITH  THE  COESa  UOLUMENI  WITHIN  1 PER  CENT  AT 
ALL  mLTITUULS  UP  TO  700  KILOMLTtHS  <2296586  FEET).  THE 
UATA  IS  AHK AuGLU  IN  THE  ATMOSPHERE  ARRAY,  A,  AS 
FOLLOWS 

A ( 1 ) = LS.  SPEED  OF  SOUND,  FT/SEC 

A ( 2 I = tl/LS)  (ULS/UZ)  , SOUND  UEH 1 V/.  1 1 WE  , 1/FT 
A(3>  = 2,  GEOMETRIC  ALTITUDE,  FT  (GIVEN ) 

a(4>  = p,  pressure.  lh/ft2 

A ( 5 ) = 0P/U2.  PRESSURE  DERIVATIVE,  LB/FT3 

A<6)  = ftHU,  OENSIlt,  SLUGS/F  T 3 

A ( 7 ) = ( 1/HHO ) (URmO/UZ I , DENSITY  DERIVATIVE.  1/FT 

A ( 6 ) = T,  TEMPERA  t URE , UEG  RANKINE 

A ( 9 ) = UT/U2,  TEMPERATURE  DERIVATIVE,  DEG  RANKINE/FT 

VARIOUS  CONSTANTS  USED 
EARTH  RADIUS  = 2089U855  FT 

SPECIFIC  MEAT  RATIO  FOR  AIR  = 1.9 

SEA  LEVEL  VALUES 

GRAVITATIONAL  ACCELERATION  = 32. 1790989  FT/SEC2 

MOLECULAR  -EIGHT  = 28.9699 

GOWMO/R*  = 0.016793916  DEG  RANK/FT 

DIMENSION  A(  9) ,HG( 10 t ,ZM(14) ,WH(19) ,TM(23) ,PM(22> 


C SET  ARRAYS  ANU  CONSTANT  VALUES 

DATA  GU,WMU,R0,GMRS/32. 1740484, 28. 9644. 20690853.0. 

1 0. 018743416/, HG/-16404. ,0.0 

2 .36UB9. ,65617. .104967. .154199. .170604., 200131., 

3 259186. ,291 ICO./, 2M/295276.. 328084. , 

4 360692. ,393701. ,992126. .524934. ,557743. ,623360. , 

5 759593, ,984252. ,1312336. ,1640420. ,1966504. . 

6 2296588. /.WM/28. 9649, 28. 88  « 28 , 56 , 

7 26. U7. 26. 92. 2b. 66, 2b. 4, 25. 85, 24, 7 ,22. 66, 19. 94. 

8 17.94.16.84.16.17/ 

DATA  I M/577, 17, 518. 6 7, 389. 97. 389. 97 ,411. 57 

1 ,487.17,487.17,454.77.325.17,325.17.379.17,469.17 

2 ,649.17,1729.17,1999.17,2179.17,2431.17.2791.17 

3 ,329b. 17. 3869. 17, 4 357. 17, 4663. 17. 4861. 17/. PM/ 

4 3711.0839,2116.2165,472.67563,114.34314, 

5 18. 128355, 2. 3162178, 1. 2321972, 3. 8U30279E-01, 

6 2. 1671 352E-02. 3. 43134 78E-03, 6. 2773411E-04, 1,53490 

7 91E> 04,3. 26242 12E*U5,1.05618G6E-05« 7. 7063076E-06, 

•i  5.8267151E-06.3,blS9854E-06.1.4S202S3E-U6.3.92905 

9 63E-0  7 ,8.40 30242E- 08.2. 2835256E- 06,7. 1875452E-09/ 
A < 3 ) = Ai 

C CALCULATE  G,  2,  ANU  CHECK 

? Z = A ( 3 ) 

G s G0*(KU/<K04Z) 1**2 

IF  (Z  ,G 1 . 295276.0)  G TO  6 

C T MS  LINEAR  WITH  GEUPOTENT 1AL,  CALCULATE  H ANU  SEARCH 

h = RO*Z/(HO+Z) 

DO  3 1 = 2«U* 

U =1-1 

IF  (HG(  I I .(it.  HI  VO  TU  4 

3 (ONTINUE 

C CALCULATE  TMS  SLOPE, IMS,  AND  SET  MOL  WT  STUFF 

4 ELM  = (TM(J+1I  - TM(O) )/(HG(U41)  - HGCJ)> 

TMS  = TM(J)  ♦ ELH* ( H - HGtu) ) 


t 


i 

■a 


1 

1 


i 


i 


i 

i 

I 


-W/J 


nt  ifGjWhiujuaiifr  I 


35 


■fggppwqg 


JTCG/AS-74-V-0 1 1 


A (2  < 1 ) S l./EM 
A12.2)  = l./TA 

A l 2 • 3 ) 2 E«U*<l./CnWA  - l./EHWF) 

A ( 2 . * ) 2 0.0 
A13.1)  2 l./EM 
A| 3.2  ) 2 l./TA 

At  3.3)2(GAMRA/CPI<HCPA-CPF-RGAS*<1./EKLA.1,/EHWF)  ) 

A13. *)  2 -CPV*TV/E 
A t * • 1 ) 2 1.0 
At*. 2)  2 0.U 
A»*.3)  2 0.0 
At*.*)  = .1.0 
Ytl)  = HUOTOG*ZOU/£H 

Y ( 2 ) = FOOT/PR*  VUOT/V 
hTHANS  2 o. 

DO  300  1=  1.3 

300  HTHANS  2 hTHANS  ♦ HJl I ) *AJl I > *(TJ< I ) - TAi 

Y 1 3 ) 2 t MOOTEV*CPF*TF  ♦ MOOTOfi*CPOG*TF  - HOOTED* tCPF*TA 
1 UtLHF ) ♦ HTHANS )/E 

Yt*J  2 HDOTEV  HOOTQG  - HUOTCU 


CALL  LLSK  <C,X.UUMnY.*.l.*.O.LA> 

IFInoarv  ,LU.  0.)  (>0  TO  37^ 

IF  t NOUTVL/MUOTV)  330  « 330  • 375 

350  NCYC  = NCYt  *1 

IF  HCYC  .LE.  10)  GO  TO  1*0 
WKX  Tt ( b . 90  0 ) 

STUP 

375  IF  (LA.EU.O)  GO  TO  *00 
WHITE  lb. 000) 

STOP 

*00  CONTINUE 

ttOO  tOHMAT  <30HCOEl  ( XCIENT  MATRIX  IS  SINGULAR) 

900  ( OHM AT (30HVLNT 1NG  CALCULATION  IS  CYCLING) 
RETURN 


mn  ono  ft  nr.  n 
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f 


75 


100 


TJ(1)  = 

L1N( 

H( 

n 

.6) 

TJ(2)  = 

LIW( 

M( 

i> 

• 3) 

TJ( 3)  = 

L1N( 

M( 

i> 

.2) 

PVAP  = 

LIN  ( 

TF 

.7 

) 

oo  loo  ; 

I = 1. 

3 

tj(i)  = 

IJU) 

* 

659 

CALCULATL  HU0T06,  HUOTCU,  HOOTEV 


HUOTCU  =0.0 
RHOKT  = PK*ENW 
PPAIH  a 2*KhQKT/EHWA 
PPFUEl  = PH  - PPAIH 
PPFL0S  = PVAP*146. 

IMPPM.8S  .GE.  U.99*PR)  GO  TO  120 
H'JQTEV  =0.0 
PPALGS  = PH  - PPFLOS 

IF  ( AbS I PP AlH -PP ALBS ) .LT  • 1 . t "10 ) GO  TO  130 
PAH  = (PPAIH  - PPALQS ) /ALOG ( PPAIR/PPAUGS) 

CKG  = UV*PH/(KGAS*1A*PAH*CUELTA) 

HUOTEV  =AJ( 3) *CKG* (PPFLBS  - PPFUEL  )*ENHF 
GO  TO  130 
120  WHITE.  (6.65) 

65  (OKHaT ( ////26H(ULL  BOILING  HOOTEV  CONSTANT) 

TtHPOHAHY  CALCULATION  OF  AMOUNT  OF  DlSOLVEOGAS 
130  EMWUG  = EMMA 

ULlg  = (GAL0-GAL >*231./172tt. 

EHUISE  = (btTA*EMWOG*VLlB  *PPAIH)/(. 797*653. *2116.220) 
EHUIS  = LMOISO  - EMOG  - SUMMOC 
IF (EMblS  .LT.  0.  ) WXITE(6.55> 

55  (OHHmT(«OAHUUNT  0(  U1SS0LVLL1  GAS  IS  NEGATIVE*) 

CON  = CONI 

IF( (tMUIS-EMDTSE)  .LT.  0.  ) CON  a C0N2 
HUOTOG  = CON.(EMOIS-EMUISE) 

EHUOOT  = EHUIS*GALUOT*231 ./VL1Q/172S. 

I'M  TIME  .LB.  TO  ) HUOTV  = MOOT 
160  IF (HUOTV  «GL • 0.0)  GO  TO  150 

VEKIlNb 


CPV  = CP 
TV  = TA 
2 V = Z 
LMUV  = LHU 
GO  TO  <fOn 


FILLING 


t 


15  -■ 


200 


C 


CPV=  CPA 

TV=  TALT  * ( 1 . ♦ SUKT(0.72)*0.2*AMACH**2) 

IF ( T VENT  .01.  0.)  TV  = TVCNT 
ZV  = 1.0 
EMWV  = EHWA 

RHu  V = PR*EMtaV/(RGAS*TV) 

MUOTVL  = MOOTV 

A (1.1)  = 2/EM 
A(l. 2)  = 0.0 
A ( 1 • 3 ) = 1.0 
A(l. 6)  = -2V/EH 
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sukjkouune  oehiv 

COHHuN/LIMITS/TO.THAA.UT.TPRNT.OTPRNI . TIME.RX6 
COMMON/STIK/  RGAS.tMdA.tMUF .LPA.CPF .HJt3) .ULL6H.ULMID.ULHT . 

XDtLHF  .ZOG.AV.CP.EMb.UV  »CU£LT  A ,K  T ANK  .bALO.TVENT 
COHMON/EULIN/  MUOT.TUOT.ZDOT.  MUOTV.  MuOTOS.  MQOTEV.  HU0TC0.W00T, 
IGALUOT  . 

2fcM . I A . Z . EH  V . EMOG , EnE  V , LMCU  . V .GAL 
COMMON/ETC'  ALT . A d ( 3 ) . Td ( 3 ) . KHOV . HR . PPAIK . PPFUEL . PPFLUS • 

X ArtACH.UALTUT.PUOT.OF .MASKAT .UVENT 
COMHON/OUTbAS/bETA .CONI ,C0N2 .KHOLlQ » LKU ISO .SUNMUO oEMOUOT 
X . VcIG.tHUIS .EMUISE 
REAL  LIN. Hill 

REAL  MUQT .MUOTV »Ml)OTtV .MUOT LU.MUOTOb.MUOTVL 
OlMENSlON  A ( A i 4 ) , Y (4 ) .CIA. 5) . X ( 4 ) .DUMMY (4 ) 

CUUI VALENCE  (A.C).  t».C(17)l.  (X.  MUOT).  ITdtSI.TF) 

NAMLLIST/Otu/  PPALGS  < Ad . T J .KHOV  « CP.DaLTUT  <CPV.TV.ZV iEHWV  <2 .ZQOT 
fJANELIST/COEF/A.  Y.hTKANS 


NCYC  s X 
M(II  = TIME 
AMACH  = L I N ( K ( X ) .8) 

CP  s (I.-2)*CPF  * Z*CPA 
EhW  s X./(  2/EHWA  ♦ U.-ZI/EMWF  i 
bAMMA  = X •/ < X • -HGAS/ (EMU*CP*77A . ) ) 
t = tM*CP*T A/GAMMA 
CPOG  s (X,  - ZOG)*CPF  4 ZOG*CPA 
ALT  s LIMHU>.5» 

ALT  s ALTaXUUO.U 
DELTs  Ut 

IM  TIME  +UT  .GE.  TMAX  ) OLLT  =-UT 
TPULL I s TIME  4UELT 
ALT X s LIM  THUELT.5  ) 

ALT!  = ALTX*XOOQ.U 
UALTUT  s ( ALT X - ALT ) /UELT 
CALL  ATMOS  (ALT ,TALT,PK«OUrtM,DUMX.DPUALT) 

PUOT  s UPOALTAUALTUT 
GALUOT  s LIN(  M(X).G) 

VOOT  s GALD0T*231«/1728. 

IM  KTANK  .EO.  X)  GO  TO  50 

KECTANGULAK  TANK 

CALC  Ad  S « TJ  S X=  TOP  . 2s  4 SIDES  , 3=  FUEL  SURFACE 

ULHT  = V/ (ULR IU*ULLGH ) 

AJ(1)  S ULb 1 0*ULLGH 
AJ ( 2 ) =2 . *ULHT* ( ULMIU  + ULLGH ) 

Ad ( 3)  = Ad ( X ) 

GO  TO  75 

CTL1NUR1CAL  TANK  - AXIS  HORIZONTAL 
X s TOP.  2 s 2 CIRCULAR  SIDES.  3 = FUEL  SURFACE,  DIAMETER  = ULUIC 

50  GALNOW  s (GALO  - GAL)  * 23X./X726, 

THLT  s «.*GALNOW/<  ULUIU**2*ULLGH) 

Add)  S (3.I4X592G5  - THET/2. ) * ULWIU*ULLGH 
Ad ( 2 ) s ULW I 0**2* THLT /4. 

Ad ( 3 ) = ULWI0*ULLbH*SIN(THET/2.) 

C 
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V 


*v 


« 


7<*  CONTINUE 
75  CONTINUE 

CALL  INI  TIL 
CALL  OER1V 
CALL  PRINT 
C 

C ENTER  INTEGRATION  LOOP 

C 

Mlfes  0 

ioo  mig=  Mi6+i 

CALL  FUNCT 

IF!  IPRT  . L8.  OJ  60  TO  20 0 
C 

C LAST  TIME  STEP  WAS  To  PRINT  STATION 

C 

150  TI«E=  1PRWI 
CALL  PRINT 

TPKNT = TPKNT  ♦ 0 TPKNT 
IF<T1ME.G£.TMAX)  GO  TO  1001 
oi=  orsv 

IPKT=0 

200  CONTINUE 

C CHECK  FOR  PRINT  STATION 

C 

CK=  TIME  ♦ OT+  TOL 

IF (T1HE.GT.TMAX)  GO  TO  1001 

IF<  CK  ,6T.  TMAX)  (PKNT=TMAX 

IF  ( CK  .LT.  TPRNT)  GO  TO  100 

IPKT=1 

OTSVs  UT 

OT=  TPKNT-  TIME 

1F(  OT  .LE.  TOL)  GO  TO  150 

GO  TO  100 

1000  STOP 

1001  GALX=GAL0 

WRITE! 7 ) (FAR ( lit  1=1 1 1255) 

GO  TO  50 
ENO 


t 
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SUBROUT  X WE,  PRINT 
COMMON/ T 1TL/CI0(12) 

COMMON/HMiTS/TO , TMAX.OT  .TPRNT  .UTPRMT.TIME.N16 
COMMON/STIK/  RGAS.EMWA.EMWF .CPA.CPF  »HJ<3) < ULLGH . ULWID « ULHT . 

10ELHF  .ZOG.AV.CP.EMW.Uv.CUELTA.KTANK.GALO.TVENT 
COMMON/EuLIN/  MOOT iTUOT «2Q0T i HUOTWt  MDOTOG.  MQOTEV « MQOTCU.VDOT. 
1GALDOT. 

2EM.TA.Z.EMV.EM0G.EHEV.EMC0.V.GAL 
COMMON/ETC/  ALT ♦ Ad< 3 ) , 1 J ( 3) .KHOV .PR .PPAIR .PPFUEl .PPFLWS, 

1  AMACh.UAE  I UT  « PUOT  .OF  , MASK AT , U VENT 
COl’iMUM/OU  f GAS /BET  A . CONI  , C0N2  . RHOLlW .EMOISO . SUMMUO .EMDOOT 
l.VLIQ.trtOIS.EnDlSE 

COMMON  /F  ILE  7/  F AK  l 250 i tFLEFT <250  > «F  TEMP  < 250 ) ,HT  <250  > , VT  <250 ) , 

11 TPKOT  , GALX  « TINC 
REAL  Mil) 

REAL  MOOT  t MOO TV  , MDUTEV  t MUOTCU « MOOT OGi MASK AT 


M<1)=  flrtE 

UVtNT  = MUOTV/<RHUV*AV) 

OF  = Z/U.-Z) 

EMI  s PK*V*EMW/<KGAS*TA) 

MASKAT=  <l.-EMl/EMi*100. 

WRIT£<fa.ll)0>  <C10<1»  *1  = 1*12) 

WKI Tt  t 6 t S00 ) M < 1 ) < PPA IK . AM ACH « PPFUEL , U VENT . PPFUJS .HASRAT . OF 
WRITE ( 6 « 6UU  > ALT  , UALTOT  »PR ,PDUT  .TA.TOOT. V * VUOT . GAL.GALDOT 
UKITE<6,700)EMV,MOUTViLMEV«MUOTEV»EMOG*MOOTOG,EMCD,MDOTCO,EM,MOOT 
lTPKNTslTPKNT+1 
FAI<<  ITPKNT  >=1.0/OF 

FLEFTUTPKNT  ) =10  0.0*  U.Q-GAL/bAEO) 

FTEMP<IIFRlJl  )=TJ<3)-459,7 

HT<ITPKNT)=ALT 

VT  < 1THKNT )=AMACH 

100  <UKMaT<1H1///33X,12A6///////) 

500  < UK  MAT  <20x,2UH  TIME  =,  6i3. 5.  27X  *23HAIK  PARTIAL  PRE 

1SSUKE  =.El3,5// 

2 2UX,20H  MACH  NUMBER  = , G13.5.27X .23HFUEL  PARTIAL  PR 

3ESSURE  = ,El3 .5// 

4 2UX.2UH  VENT  VELOCITY  =,  G13.5«27X«23HFUEL  VAPOR  PRES 

5SURE  =,E13.5// 

fa  2UX.20H1NTEGRATIUN  ERROR  =,  G13.5«27X*23H  AIR-FULL  RAT 

710  = ,E13 ,5/ 

a 20X  « 18HTOT  AL  MASS-PERCENT) 


faUO  < QPMAT  < //fe5X ,5H VALUE » 16X  , 10HUER IVATIVE// 

1 39X « 12H  ALTITUDE  , 10X , £13 . 5 • 1 OX , El  3. 5/ 

2 39X.12H  PKESSURE  . 1UX »E13 . 5, 10X.E13. 5/ 

3 39 X , 12H TEMPERATURE  . 10X.E13. 5. 10X ,E13 . 5/ 

4 39X.12H  VOLUME  » 10X »El3. 5. 10X.E13 . 5/ 

5 39 X , 12HG ALLOWS  USED , lux ,E13 . 5 «A0X *E13. 5 ) 

700  < OKMAT  < //fa4X  , 4HMASS 1 19X < 9Hl"l ASS  <LUX// 

1 40X , 1UH  VENTED  »10X.E13.5*10X»E13.5/ 

2 40X, lOHLVAPURAILU»10X«L13,5» 10X*E13»5/ 

3 40X.1UHOUTGASSEC  » 10X  .El 3.5 . 10X «E13. 5/ 

4 40 X , IOPCONOENSEU  » 10 X »£ 1 J.S . 10X .E13, 5/ 

5 40X.10H  TOTAL  . 10X .LIS .5 . 10X .E13 , 5 ) 
NAMtLiST/OUAS/LMOlSO . SUMMUO .EMUUOT « VLIO .EMUIS . EMUISE 
RETURN 
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INPUT  DATA  FOR  WSFT  B-1  FORWARD  FUSELAGE  TANK  TEST  CASE 


a 


FUSELAGE  TANK 


9 ENO  OF  kECOKO 
9 END  OF  RECORD 

TEST  CASE  FOR 

*OATA 

RGAS=1545. 

EMwA=28.9b6,EHWF=72. 

CPA=0.24,CPF=0.49, 

TA=60» 

HJ=3S2, 

ULLGH=10.0,ULl«IIU=10.0,ULHT  = G.55. 

0ELHF=1. 

Z0G=1 i 
AV=0.16, 

OV=0.3, 

CQELTA=0,01» 

KTANKsO , 

GAL0=S573. 

BETAsO.lt>. 

COlU  = 1000«CON2=0. 

TVEMT=70, 

TF (1 .1 1=0.1.2.3.4.5.6.710.9.10.11. 12 . 

TF <1. 2 I =60. 60. 56. 46. 48. 45, 15. 16 .20. 25. 110. 130. 120. 

TSIOE (1,11=0, 12, TSIOtll, 21=70, 70, 

TT OP <1.11=0. 12, TTOP <1.21=70. 70. 

AU  <1.1 1=0.. 1. .3,. 5. 1.4. 5, 6. 6. 8. 3. 12, 

ALT <1.2 1 =0.8, 20. 22, 22, 20. 18. 20, 22.. 25,. 25. 

GALOOT (1, 11=0. 8. 3, 6. 31, 10. 3, 10. 31, 12, 

GALOOT <1,2 1=0, 0.2779, 2779, 0,0. 

PtfAP <1,11=17, 41.6 7 ,96 , 129  « 166  <PVAP( 1 ,2 1 = .35« . 60.1. 1»2.0. 4. 0.8.0. 
EMI NF (1,11=0, .5,12  .EMI NF (1,21=0.. 85,. 85, 

T 0 = 0 , TMAX=11 ,0TPH!\lT  = ,25 «OT=. 00H 
9 ENO  OF  RECORD 
9 


m 
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PATH  TAPE  CREATOR  PROGRAM 


FV1.CM10000U.T20.I0100  . Pb. 

F VI. CM100000.T20. I0100, P6. 

ATTACH(TAPE7.  TNK1 ) 

REQUEST, TAPES. *PF. 

FTN(R=3) 

LGO. 

CATALOG ( TAPES. B1A.RP=999.RN=1) 

9 

PROGRAM  MAIN (INPUT, OUTPUT, TAPE5= INPUT , TAPE6 . T APE7 , TAPE8 . TAPE9 , 

IT  APE 10  « T APE1 1 »T  APE12  .TAPE13.TAPE14 « TAPE1S « TAPE16 ) 

COMMON  /SAVE/  HT<250 ) ,VT (250 ) ,FAR(250.10 ) ,FLEFT(250 *10)  i 
IF  TEMP (250, 10 ) .GAL (10) .NT ANKS .NTPO , T INC 
COMMON  /FILE7/  FA(250 ) .FL ( 250 ) .FT (250 > .H (250 ) . V (2S0 l .NTP.GALX.TIN 
INTEGER  GAL 
K£AO<5, 1000) N TANKS 
REAO(7)(FA(I). 1=1. 1253) 

00  1 1=1, NTP 
HT ( I ) =H ( I ) 

VT ( I ) =V ( i ) 

FAPti.l)=FA(l) 

F TEMP ( 1 « 1 ) =F  T(  i ) 

1 FLEF  T (I.1)=FL[1) 

NT=NTP 

NTPOsNTP 

T1NC=TIN 

OAL ( 1 ) =G ALX 

IF  (NTANKS.EU.l)  GO  TO  4 
1EN0=6+N TANKS 
00  3 1=6.1ENU 
K = I-6 

READ (I ) (FA«d) .J=1.1253) 

IF(MT.NE.NTP)  GO  TO  5 
00  2 J=1 « NTP 

1F(M1 (d) .NE.H(J) ) GO  TO  5 
IF ( VT (d)  .NE.V( J) ) GO  TO  5 
FAR(d.K )=FA(d) 

FTEMP(d,K)=FT(d) 

2 FLEFT(d,K)=FL(d) 

3 GAL ( K ) =GALX 

4 WftlTE(G) (HI (I) . 1=1.8013) 

STOP 

5 PRINT  2000 
STOP 

1000  (ORMAT(IS) 

2000  (ORM*T<*  01 ( (ERENT  (LIGHT  PROIILES  (OR  TWO  TANKS*) 

END 

9 

1 

9 

9 
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VULNERABILITY  PROGRAM 


• » 


i 


t 


i 


» 


> 


ms* 

m si,.,,. 


FVM  .CM60000  «T600  < 10100 ,P2. 

FTN(R=3,H=FVMAO) 

ATTACH ( TAPE6.F9U) 

FWMAD. 

9 

PROGRAM  F VHN  ( INPUT « OUTPUT  .Tf.PE.fa) 

COMMON  /COUNTER/  VSCUM ( 1U ) , RSCUM (10  > 1 1HI T < 10 ) . IEF (10 ) , XAFEN ( 10 ) * 

II AFEXA( 10 ) »IAFEXb(10) , XRAM (10)iIFE(10), I LEAK ( 10 > 

COMMON  /WEAPONS/  VZERO ( 32) ,PENM( 32) t CORAG (32 ) » ISORF (32) , AZHIN(32) • 
1AZMAX (32),ELMIN(32) * ELM AX ( 32 ) .NSHOTS ( 32 > , TIGN( 32 ) . TASH ( 32 ) « IHE( 32 ) 
COMMON  /TANKS/  DSKIN(6«10 ) .1X6(6,10) «XCG (10 ) , YC6( 10 ) . ZCG( 10 ) « 
1X1U0) ,Y1(10> ,Z1(10) ,X2(l0>  ,Y2(10>  ,Z2(lO>  ,H«IN(10>  .FULLPCUO) 
common  Pinout/  a,b,c, i asp, jasp.xin, yin, zin.xqut, tout, zout, entry, 
iafen,afex,xo,yo.zo 

COMMON  /PROFILE/  H(250) « AMACH(2S0) ,FaR(250«10) « PFUEL ( 250 < 10 ) « 
1FTEMP (250 , 10 ) tGAL(10) »NTANKS»NTP»TINC 
COMMON  /ALPHA/  AT ANK ( 2 . 10 ) . AWEAP ( 32 > . AWAR ( 2 ) , ASOL ( 2 > , ATARG 
COMMON  /M1LERRS/  EMIL(32,32) 

COMMON  /FRAGS/  VE(32) ,PSX ( 32 ) .COSMAX ( 32 ) .COSMIN (32) ,CD(32 » ,$F (32 » « 
1XF(32> ,YF(32) ,ZF(32) 

DIMENSION  CEF(IO) .CFE(IO) ,CLEAK(10> .PHCUMI 10 ) .RAF ( 10> ,PBF(10) 
EQUIVALENCE (CEF (l)«XEF(l))«(CFE(lt,IFE(l)), (CLEAK ( 1 ) « 1LEAK (1 ) I « 
KPHCUM(l)  tlAFEXA(l) > 

LOGICAL  ENTRY. NOHE.AFEN.AFEX 
INTEGER  GAL.XST«W£APS(250) 

DATA  TWOP1.UTOR, GEE/6. 26318531, 0,017953293. 32. 172/ 

DATA  AWAR , ASOL/IOHFRAGMENTIN. 10HG  WARHEAD  .10H  SOLIO  SHO.IOHT  REAP 
ION  / 

CALL  REPT1 

READ  20Q0,NwEAPS,XST,HKAM,KATMIN,RATMAX 
READ  2001 . ( WE APS (II. 1=1 .NTP ) 

REAU  2002 . ( ISOKF ( IT , 1HE ( 1 ) .NSHOTS ( 1 1 , AZMIN ( I ) . AZMAXII 1 .ELMIN ( I ) , 
1ELMAX (I).XF(I) ,YF (I) ,ZF(I) ,I=1.NWEAPS) 

REAU  2003 . ( T IGN ( I ) , T ASH ( I ) , VZERO ( I ) ,PEnM ( 1 ) .CDRAG ( 1 ) , 1=1 .NWEAPS ) 
READ  2007,(VE(i).PSX(I) » COSMAX ( X ) « COSMIN ( I )*CO(I).SF(I)t 
11=1, NWEAPS) 

READ  2009, ( (EMIL (I, J) »I=1,32),J=1 .NWEAPS) 

READ  2005, (UIG(I,d). 1 = 1, 61»(0SKIN(I,U). 1=1, 6>»J*1.NTANKS> 

READ  2006,  (XCGU)  ,YCG(  I ) ,ZC6(  I ) ,HMIN(1)  .FULLPCd)  ,I  = 1,NTANKS) 

READ  2007,(X1(I),Y1(I),Z1(I) ,X2 ( I ) , Y2 ( l ) ,Z2 ( X) » I«1 ,NTANKS) 

READ  2006, (AWEAP(I) ,X=1, NWEAPS) 

DO  103  ITIMEsl.NTP 
TIM=TXNC*(XTIME-1) 

VT=0EMACH(AMACH(ITIM£) ,H(ITIHE) ) 

CELT  A=0 , 5/VT 
ALT=H(XT1ME) 

DO  102  IWEAPsl, NWEAPS 

IF ( AND (WEAPS(ITIME),2*#(IWEAP»1)),NE,2**( I WEAP >1 ) ) GO  TO  102 
SLUGS&PENM ( I WEAP ) / ( 7000 , U9GEE ) 

VEMsVE ( lwEAP ) 
b6=CU(IWF.AP) 

COSTMX=CO$MAX ( XWEAP ) 

COSTMN=COSMIN ( IWEAP ) 

PSIST=PSX ( 1 WEAP ) 

FMsPENM(IwEAP) 

SIGF=SF( IWEAP) 

XFUSE=XF( IWEAP) 

YFUSE=YF( IWEAP) 

ZFUSE=ZF ( IWEAP) 

T1=T1GN ( IWEAP ) 
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T2=TASH< XWLAP) 

XVlsVT 
I=IVT/100+1 
XFU.GT.3Q)  GO  TO  . 

tKRbtMxm  ,iweap>  + io.oi*vt-i+i.o>*(ewil(i+i.iweap>-emil<i,iweap> ) 

00  TO  2 

1 tKKstMXU X *1WEAP) 

2 NOHL=.TRUE. 

IF  UHLlXwtAP)  .OT.O)  NOH£=. FALSE. 

NSiHsNSHOT S ( IWEAP  ) 

00  500  1=1,100 
500  VSCUMU)=0«0 

00  100  lSHOT=l,NSH 

A2=AZMIN ( IWEAP ) +R  ANF ( XST ) * ( A2MAX ( I WEAP ) -AZHIN ( IWEAP ) ) 

£L=ELMXN( IwEAP)+RANF (XST  )*  (ELMAX  { IWEAP  > »ELMN(  IWEAP ) ) 

AZ=AZ*DTOR 

EL=EL*DTOR 

COSE=COS(£L) 

XMU=-C0SE*C0S<AZ> 

VHU=-COSE*bXN(AZ) 

ZMU=S1N<EL) 

KS=AbS( ALT/ZMU) 

XH ISORF ( XWLAP) .GT.l ) GO  TO  3 
SIGMA=KS*EKK*0.0  01/0.£>745 
GO  TO  4 

3 SIGMA=ERK 

4 0M=A13S  (GAUS(SIGMA)  I 
THETAsHANF (XST)*TWOPI 
COSTH=COS(THLTA) 

31i\lTH=SIN(  THETA) 

VSHsVZLKO ( XWLAP ) *EXP ( -RS*CORAG ( IWEAP ) J 
AA=VLIVVSH 

VK=SOKT(VSH*VSH+VT*VT-2.0*VSH*VT*XMUJ 

As(OSH*XMO-VT)/VR 

B=VSH*YMU/VK 

C=l.U-A*A~b*b 

IF ( C ) 5,6 « 7 

5 B=SIGN{5URT(1,0-A*A) ,YMU) 

6 C=0,U 
GO  TO  a 

7 C=SIGN(SWRT(C),ZMU) 

IFtABS(C) , to. 1.0)  GO  TO  9 
UsSQKT ( A*A+b*B ) 

COSPSX=-A/U 

S lNPSl=B/0 

XO=OM*(C*COSTH*COSPSI+SINTH*SINPSX) 
YO=OM*(-C*COSTH*SINPSI+SINTH*C0SPSI ) 

Zu=(-A*XO-b*YO)/C 

go  ro  10 
a X0=DM*B*SINTH 
Y0=-0M*A*SINTH 
ZO=OK*COSTH 
GO  TO  10 
9 XO=OM*COSTH 
YO=UW*SINTH 
Z0=0.0 

10  IF ( ISOKF (IWEAP ),GT«1)  GO  TO  50 
ENTRY=. FALSE, 

00  11  ITANK=1,N TANKS 

CALL  IN( XCG(  I TANK ) » YCG ( I TANK ) ,ZCG( I TANK ) ,X1 ( ITANK ) * Y1  ( ITANK  ) « 
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1ZK  I T ANK  ) , 0 . 0 1*PFU£L  ( I T IME 1 1 TANK ) *FULLPC  ( ITANK ) ) 

1F( .NOI .ENTRY)  GO  TO  11 

CALL  OUT  <XCG( ITANK) » YCS< ITANK) ,ZC6< ITANK) *X1 ( ITANK) «Y1( ITANK) 
1 Z 1 ( I T ANK ) , U . 0 1 *PFUtL ( I T I ME « ITANK) *FULLPC (ITANK)) 

GO  TO  14 

11  CONTINUE 

IF (NONE)  GO  TO  100 
00  12  ITANK=1,NTANKS 

CALL  lN(XCGUTANK)  «YCG(  ITANK ) »ZC6( ITANK ) .X2C  ITANK ) «T2C ITANK )» 
1Z2CITANK) .-1.0) 

IF(. NOT. ENTRY)  GO  TO  12 
ILEAK ( ITANK ) sILEAK ( ITANK )♦! 

IHIT(ITANK1=IHIT(ITANK)+1 

IF(I1G(IA$P.ITANK).6T,0)  IEF ( ITANK )=IEF ( I TANK) 41 

12  CONTINUE 
GO  TO  100 

14  IHIT ( ITANK )=IHIT< ITANK)+1 
RATIO=FAR(ITIME. ITANK) 

RSCUM(ITAi|K)sRSCUM(ITANK)+RS 
VSCUM  ( ITANK  )sVSCUM( ITANK )+VSH 

ENEKGY=0.5*^LUGS*VSH*WbH 

01=DSKlN ( I ASP .ITANK ) 

IF ( ENERGY «LT ,HM IN (ITANK) ) GO  TO  100 

IF(AFEN)  GO  TO  23 

IF (ENEKGY.LT .HRAN ) GO  TO  15 

IRAK < ITANK )=IRAM( I TANK) +1 

GO  TO  100 

15  1F(AFEX)  GO  TO  17 

ILEAK (ITANK) sILEAK ( ITANK )+l 

IFUIGIJASP. ITANK). GT.O)  GO  TO  16 
200  IFdiGUASP. ITANK). GT.O)  GO  TO  16 
IF ( VSH*T 1.GT.01)  GO  TO  100 
IF(VSH*T2,L1 ,01)  GO  TO  100 

16  IEF(1TANK)=1EF(1TANK)+1 
GO  TO  100 

17  IAFEXB(ITANK)=IAFEXB(ITANK)+1 
IF (RATIO. Lf.RATHIN)  GO  TO  22 
IF (RATIO. GT.RATNAX)  GO  TO  22 
GO  TO  (16, 18, 19. 19.20.20) *1ASP 

1#  02=01*2. 0*X1 (ITANK) 

GO  TU  21 

19  02=01+2. 0*Y1(ITANK) 

GO  TO  21  4 

20  02=01+2. 0*Z1 ( ITANK ) 

21  IF(VSH*T1.GT,02)  GO  TO  22 
IF ( VSH*T2.LT ,01 ) GO  TO  22 
1FEUTANK)  = 1FE  (ITANK  » + l 
GO  TO  100 

22  ILEAK (ITANK ) sILEAK (ITANK ) +1 
GO  TO  200 

23  IAFEN( ITANK )=IAFEN(ITANK)+1 
IF (RAT10.L1 .RATMIN)  GO  TO  26 
IF (RATIO. Gi .HATNAX)  GO  10  26 
GO  TO  (24,24,25,25,26,26) tIASP 

24  02=01+2. 0*XI(ITANK) 

GO  TO  27 

25  02=01+2. 0*Y1 < ITANK ) 

GO  TO  27 

26  02=01+2,0*21 (ITANK ) 

27  IF  <VS>H*T1.G!  .02)  GO  TO  26 
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IF (VSH*T2,LT,01)  GO  TO  26 

1FE  ( 1 TANK ) = 1F  E(  HANK  ) +1 

IF  ( AFEX)  IAFLXA(ITANK)=IAFEXA(ITANK)*1 

GO  TU  100 

26  IF  ( .uOT  • AFtX  | GO  TO  29 

I AF  EX A ( 1 TANK ) =1 AFEXA ( I TANK ) +1 
GO  TO  100 

29  IF (ENERGY, LT.HKAh)  GO  TO  30 
IRAKI  1 f AUK ) = IRAM( I TAWKJ+l 
GO  TO  100 

30  ILL AK ( I TANK  ) =ILEAK ( I TANK  ) 4-1 

lFtllGlJASiF,  I TANK  ) . GT.Q ) 1EF { 1TANK > =IEF ( 1TANK > +1 
GO  TO  100 

50  TF=GAUS(SIGF »/VSH 

XST  AR=XO+TF* ( VSH*XKU-VT ) *XFUSE 
YSTARsYO+TF  *VSH*YKU+YFUSE 
ZSTAKsZO+TF*VSH*ZMU+ZFUZL 
00  61  ITANKsi.NTANKS 
YT  = YCG( ITANK  J-YSTAR 
ZT=ZCG< IT  AUK ) -ZSTAR 
XTT=XCG( I TANK l-XSTAR 
C1=VR*WR-V£K*VEM 

C2=0SH* ( (XMU-VT /WSM) *XTT  +YMU*YT ♦ZKU*ZT ) /Cl 

C3-X|T*XIT+YI*YT>ZT*ZT 

()ISCR=SuRT(C2*C2-C3/Cl) 

T=C2-UISCR 

IF(T.LT.U.U)  T=C2+01SCr< 

ITER=0 

51  IF (ITLR.6T.25)  GO  fO  61 
XT=*7T-*0T*T 

XL=SURT (XT*XT*YT*YT*ZT*ZT) 

IF(XL.EQ.O.O)  GO  TO  61 
BLT  AX  = XT/XL 
BETAY=YT/XL 
BLT  AZ=ZT / XL 

COSGAM=XMU*6LTAX+Ymj*BETAY+ZMU*BCTA2 

VO=Vbh*COSGArt+SURT (USH*YSH*COSGAM*COSGAM+VEM*WEW-VSK*VSH) 
VOOOT=VO*VSH*VT*(XMU-BETAX*COSGAM)/(XL*(VO-VSH*COSGAM>) 
FOOT=BO*( ( 00+VODOT*T)/(l,Q+Bb*V0*T)-BETAX*VT> 

IF (FOOT. LE. 0.0)  GO  TO  61 

F=AL0G(1,Q+BD*V0*T)-6B*XL 

TNEW=T-F/FDOT 

IF (ABS(TNEW-T) .LE. DELTA)  GO  TO  52 

T=TNEW 

ITEftslTEFUl 

GO  TO  51 

52  IF ( XL .GT . 500 ) GO  TO  61 
COSTH=( VO*tOSGAK-VSH(/VEM 
IF(COSTH.LT.tOSTKX)  GO  TO  61 
IF  (COSIH.GI  .COMMN)  GO  TO  61 
IHIT(ITA1mK)  = IH1T(ITARK)+1 

E=AA*AA*ABS(AA+COSTH)/(AA*AA  + 2. 0*AA*COSTh+l . 0 ) **1 . 5 
PS1C=PS1ST/(E*XL*XL) 

VHI  T=VO*EXPt-BB*XL) 

VNET=J>QRT  ( VH1T*VHIT-2.0*VT*VHIT*8ETAX*VT*VT) 

EOENJjsO  .5*SLUGi>»Vf4ET *VNET *PSIC 
IF (EOENG.GE.HRAM)  IRAK ( IT ANK J = IRAM ( ITANK ) +1 
ATB  = 4.0*X1UTANK)*Y1(ITANK) 

ASS=4.0*A1 (ITANK )*Z1  (ITANK) 

AFR=4.0*Y1 ( IT ANK)*Z1( ITANK) 
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SARsAUS ( 8ETAX  ) * AFR+ABS (BET AY ) *ASS 
TAR=SAR-fABS(BETAZ)*ATB 

PHCUM  < X TANK  > sPHCUM ( I TANK ) +1 . O-EXP  < -PSI C*T AR  ) 

ARBF  =0 . (U*PF UEL ( 1 T iMfc , I T ANK ) *0 . 0 l*FULLPC ( X T ANK ) *SAH 

ARAF=SAR-AR8F 

IF(8ETAZ»55.55.54 

53  ARAF=ARAF-8ETAZ*ATB 
60  TO  55 

54  ARBF=ARBF>BeTAZ*ATe 

55  PHITBF=1.0-EXP(-PS1C*AABF ) 

PHXTAf  =1. 0-EXP  (-PSIOAKAF) 

TMPsFTEMP ( IT1ME* ITANK ) 

IFtALT.LI .10000.0)  SO  TO  Sfa 
AF  ACs2 .5*EXP  ( - , 000()92*ALT ) 

IF (ALT. ST .60000. U ) bO  TO  SB 
1FOMP.LT. U.O)  60  TO  58 
IF(TNP.GT.45.0)  bO  TO  57 

DF  = TMP/45.Q*U.2-0.0Q002*ALT) 

60  TO  bU 

56  AFACsl.O 

IFOKP.LT. 0,0)  60  TO  58 
IF (TMP.GT .45.0)  60  TO  59 
DF=TMP/4S.O 
GO  Tu  6U 

57  UF=1.2-0.0UU02*ALT 
GO  TO  6U 

58  OF  so  , 0 
GO  TO  60 

59  OFsl.O 

60  CEF ( I TANK ) =CEF ( I T ANK ) +0 . 3*QF*PHIT8F 

CLEAK ( ITANK  >=CL£AK ( ITANK ) +PHI TBF* ( 1 , 0-0 . 3*DF ) 

PAT  IO=FAK(  IT  IME.,  ITANK) 

IF (KAT10.LI .KA1MIN)  bO  TO  b2 
IF (RATIO, bl • K A t MAX  I GO  TO  62 

CFt < I TANK ) =CFE < IT ANK ) ♦ . U0000769*S8KT t FM> *VNET*AFAC*PHITAF 
62  VSCUK ( ITANK ) =VSCUM ( IT  ANK ) +VNET 
KSCUMU  TANK  )=RSCUH(  ITANK  )+HS 

61  CONTINUE 
100  CONTINUE 

IF (ISOKF ( IWEAP)  ,6T.l)  GO  TO  108 
UO  101  ITANKsl.NTANKb 
IF(IHXTdTANK)  ,EU.U)  GO  TO  101 
VSAVt=VSCUM{ ITANK )/IHIT( ITANK) 

PAFEnls  (1. 0*1  AFLN(  ITANK)  ) /IHIT { ITANK ) 

IF ( IAFEN( ITANK) .E8.0)  GO  TO  104 
PAFEXA=(1.0*IAFEXA( ITANK) J/IAFEN (ITANK) 

GO  TO  105 

104  PAFEXA=O.Q 

PAFEXB= (1. 0*1 AFEXB< ITANK )) /IHIT (ITANK) 

GO  TO  107 

105  IF(IAFEN(ITANK)»EQ.1HIT(ITANK))  GO  TO  106 

PAF E XB= ( 1. 0*1 AFEXB < I TANK ))/( IHIT (ITANK > -I AFEN( ITANK) ) 

GO  TO  107 

106  PAFEXB=Q.O 

107  PLEAK=(1.0*ILEAK(ITANK)-1EF( ITANK) ) /IHIT ( ITANK ) 

PFE=( 1.0*IFE( ITANK ) )/IHlT( ITANK) 

PEF=(1.0*IEF (ITANK) ) /IHIT (IT ANK) 

RSAVL=RSCUM( ITANK) /IHIT (ITANK) 
phit=(I.o*ih;t(ITank) )/nsh 

PRAK=(1.0*IRAM(ITANK) )/ IHI T ( ITANK ) 
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(•'NueF  =1  .O-PLE  AK-PRAM-PFE-PEF 
PbFEN=1.0-PAFEN 


P8F£XA=1.0-PAFtXA 

PBFEXU=1.0-PAFtXb 

PRINT  1000*  ASOL.ATAHG.  ( AT ANK  ( I , X T ANK  ) * I = 1 * 2 > 

PRINT  lOul.TlR.AWEAP(IWEAP)  .PBFEN 
PRINT  1 002. PHI 1 »PFUEL( lTlMt.il ANK) 

PR1N1  1003.VSAVt,FTEMP( l TIME, IT ANK) , PbFEXB 
PKI.MT  10JH.KSAVE.FAR1  ITImE.ITAUK)  . PAP  EXB 
PH  111  I 1 0 0 b , H ( I TIME  ) . V T . HaFLN 
PK1NI  100t>»PBF£XA.PAFEXA 
PRIM  1007.PNOEF  , PLE  AK  . PtF . PHAM  » PFE 

101  CONTINUE 
00  To  1 02 

lUa  DO  109  I fANK=l .NTANKS 

IMPhCOMlIl  ANK)  .to. 0.0)  GO  TO  100 
VbAVE=VStUM(  1TANk)/IHIT( 1TANK) 

PLEAK=CLLAK( I TANK 1/IhITl 1TANK) 

PFE=CFE( 1TANK )/IHIT  ( IT  AUK  ) 

PEFsCtF  ( ITANM/lHITl  IT  ANK) 

RSA VEsRSCuM ( 1TAMK1/1HIT  1 1 TANK ) 

PHI T=PhCUH< IT ANK ) /NSH 

PH AHs (i.O* IK AMI  I TANK ) l/IHIT ( I TANK  1 

PRIM  1000. AwAK.ATAKO, ( AT ANK ( 1 , 1 T ANK ) .I=lt2> 

PRINT  1 0 1 3 . T 1M . AWLAP ( I WE AP ) 

PRINT  1008, PHIT .PFOELIITIME. IT  ANK ) 

PR  I Ml  1009.VSAVE.FTtMP( ITIME.ITANK) 

PRINT  1010 . SIGMA, F AH (1TIKE.1 TANK) 

PRINT  lull. ALT. VT 

PRINT  1012. PLEAK.PEF .PRAM, PFE 

109  CONTINUE 

102  CUNUNUE 

103  CONTINUE 
STOP 

1000  lOKMAT (*1*.51X.*«UEL  TANK  VULNERABILITY  M0DEL*/62X . *REPOHT  2*/*Q*. 
155X ,2 A10/*0*. H9X,* VEHICLE--*, A10/*0* » H9X , *FUEL  TANK— *.2A10> 

1001  (OR«AT(*OT1ME  INTO  MlSSION(HKS) —*. ( 7 . 3 » X8X . *TH8E AT  — * . A10 . 21X . 
l*PKOUABILITY  OF  LIUUID  ENTR  Y-*  ♦ F9 . 6 ) 

1002  (OKrtATUQPKObAblLIlY  01  HIT  ON  (UEL  TANK — *,(9.6,7X,*PERCENT  (UEL 
lREMAl NlNG--* , F7 .2  »8X , *GI VEN  LIUUID  ENTRY--*) 

1003  (OK'jIaT(*OAVERAGE  STHIKINS  VELuCITYl (PS)--*, <B,1,9X.*(UEL  TLMPERATU 
IKE ( F ) --* , F 7 «2 . 15X » *PROBAblLl TY  OF  LIUUID  EXIT— *,F9,6  ) 

1004  ( UKHA  T 1 *0AVEKAGE  SLANT  RANGE  ( ( T ) — * . 1 9 . 1 » 15X , * ( UEL/AIR  RATIO—*. 
IF9.G,18X,*PKO0ABILIIY  OF  VAPOR  EXIT— *,F9.6) 

1005  (OKVIAT<*OAIKCRA(  T ALTITUDE  1 < T ) --*. < 8. 1/*0AIKCRA ( T SPEED! IPS)--*. 
lF7.I,60X.*PR0BAbILITY  OF  VAPOR  ENTRY  — * .F9.6 ) 

lOUfc  ( OP  MAT ( *0* , 89X , *61 VElv  VAPOR  ENTRY  — */*0*,93X, *PROB ABILITY  01  LI8UI 
ID  EXIT— *,F9.6/*0*«93X,*PROBAaiLITY  OF  VAPOR  EXIT— *,F9.6) 

1007  <OPNpTl*-*,49X»*PRObAHlLlTiES  0(  (UEL  TANK  DAMAGE*/*0*.60X»*&IVEN 

1A  HlT*/*0*«5lX,*P(NO  EFFECT)  =* , F9.6/*0* , 51X ,*P ( LEAK 

2 N 1 1 HOUT  FIRE)  =*,F9.6/*Q*,51X.*P(LEAK  AND  EXTERNAL  PlRE ) r*. 

3F9.6/*Q*»51X,*P(UEsrK0CTIVE  RAM)  =*,F9.6/*0*,51X, *P( INTERNA 

4L  F IKE /EXPLOSION ) s*  «F  9 .6 ) 

1008  (ORNAT l *OPROBAbILIT Y 01  hIT  ON  (UEL  TANK--*. (9.6. 7X,*PERCENT  (UEL 
1REKA1NIU&--*,F7.2) 

1009  (OPNA1 (*0AVE|(AGE  S ( R IK ING  VELOC IT Y ( ( PS  ) — * . ( 8 . 1 , 9X , * ( UEL  TEMPERATU 
lKt(F)—  *,F7.2> 


1010 

1011 


(ORIIAT (*0AIMING  SIGrtA < ( T ) --* , < 8 . 1 , 25X . * ( UEL/ AIR  RATIO--* ,(9.6) 

( OHM AT l*UAIRCRA(T  ALTITUDE ( (T)--*» ( 8 . I/*0 AIRCRA ( T SPEED! (PS)—*. 


IF  7. 1 ) 
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1012  ( ORRAT  ( » 0*  >49X  ,*PR0BAB1LIT1ES  0(  (UEL  TANK  UAMA6E*/*0*  »60X.  *GI  VEN 

1A  H1T*/*0*«51X,*P( LEAK  WITHOUT  F IRE  ) s* ,F9.6/*0* .51* . *P ( LEAK 

2 UNO  EXTERNAL  FIRE ) = , , F 9. 6/*U* , 51X . *P ( UESTKUCTI WE  KAMI  =*, 

3F9.6/*0*  .5IX, *P< INTERNAL  F IKE /EXPLOSION ) 2*  « F9 , 6 ) 

1013  <OR**H*071ME  INTO  MISSIONIHRS I --*•  < T. 3 • 1BX  .-THREAT A10 ) 

2*J0u  ( OP  -'A  T ( 2 1 1 0 . 3 ( 1 0 . U ) 

20ul  (URrtAf) lOIfl) 

2002  (0P.1AU2U«I«i  7(10.0) 

200  5 (URiAf  (3(10.0). 

200*1  ( OR‘1  A T (8(10.0) 

2 005  I0»*6t (4X, 611. 6(10. 0) 

2006  < ORB  A T ( 5 ( 1 0 • 0 ) 

2007  <Ofti'iAT<6<10.0) 

20 UH  (0 RAAT(HAIU) 

EmO 


SUBROUTINE  IN (XT. T f • ZT , LT , WT , HT « PFUEL ) 

LOGICAL  ENTRY. AFEM 
REAL  LT 

CUnnoN  /1NOUT / A.B.C.IASP, UASP.X1N. YIN, 21N.XOUT.YOUT, ZOUT, ENTRY, 
1AFEn.AFEX.XO, YU, 20 
1 F ( C ) 1 , 4 . 2 

1 1ASP=5 
ZlNsZ  f +HT 
GU  TO  3 

2 IA5P=6 
2IUSZT-MT 

3 YINSb/C* ( ZIN-20  ) irO 
XIN=A/C*»  (2IN-201+X0 

IF(  A0S(  Ylr.-YT  ) .GT.WT)  GO  TO  4 
IF«ABS(X1N-XT) .GT.LT)  GO  TO  4 
GO  TO  12 

4 IF<0)5.O.6 

5 IASP=3 

T IN=Y  T +KT 
GO  TO  7 

6 IASP=4 
YinbYT“WT 

7 2IM=t/B*(YlN-YUU20 
XlN=A/B*(YlN-YO)+XO 
1F(ABS<Z1N-ZD.GT.HT)  GO  TO  S 
IF(AQSCX1|,-XT). GT.LT)  GO  TO  8 
GO  TO  12 

6 IF (AI9, 17.10 
9 IASP=1 
XIN=XT+LT 
GO  TO  11 

10  IASP=2 
XINsXT-LT 

11  ZIN2C/A*(X1N-X0)+Z0 
YIN=B/A*(X1N-XU)+Y0 

IF (A8S<21N-ZT) .GT.HT)  RETURN 
TFlAfciS(TlN-Yn.faT.GT)  RETURN 

12  ENfRYs.TRuE. 

1 F ( PFUEL )l7, 16.13 

13  XMPFUEl.EQ.100.0)  GO  TO  15 
1F(IASP»5)14, 16.15 

14  IF (50. 0*CZIN-2T+HT)/HT.GT. PFUEL)  GO  TO  16 

15  AFEn=. FALSE. 

RETURN 

16  AFENs . TRUE • 

17  Return 
End 
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SUBROUTINE  UUT(Xr,YT.2T,LT.WT.HT  «FF  UEL  I 
LU&ICAL  AFEX 
PEAL  LT 

CO«rON  /I MOOT/  A.B.C.IASP.JASP.X1N.Y1N-  UN » XUUT , YOUT , ZOUT .ENTRY , 
lAfEN.AFEX.XO, YO.ZO 
IF  TOl.4.2 
1 jaSP=G 
ZOUT=ZT-MT 
GO  TO  3 
i JASP=5 
ZuU  T =ZT ♦HT 

i YOUT=ts/C*tZOUT-ZO)FYO 
XUUI=A/C*UOUT-ZU)+XU 
IF l ABS1 YOuT-TT ) ,gT .WT ) GO  70  4 
If l ABS TXOUI-XD.Ol.tl)  GO  TO  4 
GO  To  13 

4 IFtblb.h.b 

5 JASP=4 
YOU  T = Y T -AT 
GO  TO  7 

G JASP=3 
YOUTsTT+HT 

7 Z0UT=C/B*1 YOUT-YOI+ZO 
X0UT=A/H*1 YOUT-YO )+XU 
IF TABSTZOUT-ZT | .GT.HT ) GO  TO  0 
lFTA8STX0Ul-xD.Gr. LI)  GO  TO  6 
GU  TO  13 
4 IF TAJO, 12,10 
9 JASP=2 

XOUT=XT-LT 
GO  TO  11 

10  oASPsl 
XOUT=XT+LT 

11  ZOUT=C/A*TX(JUT-XO)+ZO 
YOUT=b/A*lXOUT-XO)*YO 
IFTAQSTZOUT-ZT) .GT.HT)  GO  TO  12 
IFTAeSTTOUI-YD.GT.WD  GO  TO  12 
GO  TO  13 

12  PRINT  1000 
STOP 

13  IFTPFUEL.EQ.0.0)  GO  TO  15 
IF  (PFuEL.EO. 100.0)  GO  TO  14 
IF  UASH.EQ.fa)  GO  TO  14 

1FT5U.O*TZOUT-ZT+HT l/HT.GT.PFUEL)  GO  TO  15 

14  AFEX=. FALSE. 

return 

15  AFE  X = , TRUE . 

return 

1000  TUR'IATT*  entry  cut  no  exit*) 
end 

F UNO T I ON  UEFACHT AMACH.HT) 

IF THT.LT .3&069.0  1 GO  TO  1 
OE*Al.Hs9fa8.452*AWAtH 

Return 

1 UEi**  ACHs49  , U4u  772*SttRT  (51f*. 686-0. 0035661fa*HT)  * ATTACH 
Return 
END 
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SUBROUTINE  REPT1 

COMMON  /PROFILE/  HI (250 ).VT(2bO) .FAR (250.10)  .FEEFT  (250  * 10 ) « 
IF  TE«P ( 250  «10t  <GAL  ( 10)  .NTANKS.NTP.  T INC 
COMMON  /ALPHA/  AT AHK (2 • I(l ) • AHEAP (3a) . AWAK  <2 ) . ASOL (2 } .ATARG 
INTEGER  GAL 

RE At  (6)  <HH  I)  .I  = l.a013) 

RE AO  2000  «( ( A TANK (I.U) .1=1.2) < J=1 .NTANKS  > .ATARG 

NLEFTsNTaIiKS 

1ST  =1 

1 IENOslST+4 
NLEF  T -WLEFT  «5 

IFi. JLEF  T. 61 .0)  GO  TO  2 
XEiVJslEIJU+uLEFT 

2 LINE=0 
t',TP.i=l 

2 IF  I MUU (LINE  . GO ) .GT • 0 ) GO  TO  A 

PRINT  1000.  ATARG,  ( (ATANMd.I)  .J=1.2)  .I=IST.IENU) 

PH  I **'l  T 1001.(GAL(1)»I=IST«I£NU) 

PRINT  1002 
L IN£=LINE+12 
4 Tir=lluc*(NTPS-l» 

PRIM  100, 3* TIM. (FAR(UTPS.I) .FLEFT (NTPS. I ) «I=1ST.IEND) 

L 1NE=L JNC+1 


MTPS=,jTPB*l 

lFTi'JTPS.LE.NlP)  GO  TO  A 
IKNLEf T.LE.O)  RETURN 
1ST  =1ENU+1 


(.0  TO  1 

1000  ( OK HaT (*1*»51X«*( UEL  T AUK  VULNERABILITY  M00EL*/62X • •REPORT  1*// 
15»X  .* VLhitLE—  *«A10/1  IX. ST4X.2A10)  ) 

1001  (URnAT(llX.5<SX.2H*-*)»/llX.5(l8.*  GALLON 


1002  (UR.'IAT  (liX.5(3X.21(*-*H/*  TIME 


1*  MISSION  • « 5 ( * 

1003  < OKMATT IX. T 10. 5. 5 ( (12.8.(12.2)) 
2000  (UHMAU5A10) 

ENO 


F/A  RATIO  REMAINING*) ) 


CAPACITY*)) 
INTO* • 5 ( 15X . *PCT • (UEL*)/ 


FUNCTION  GAUS ( SIGMA ) 
CAuS=0.0 
00  1 1=1.12 
1 GAUS=GAUS+RANF(0.0) 
GAUS=SlGi1A*  ( GAUS  -6  • 0 I 
RETURN 
END 

9 END  OF  RECORD 


51 


ffryg rrsgggg  yaesrzs  aaaasa— te- 


JTCG/  AS-74-V-0 1 1 
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Attn:  SAVDL-EU-MOS  (J.D.  Ladd) 

Attn:  SAVDL-EU-MOS  (C.M.  Pedriani) 
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Army  Aviation  Systems  Command 

P.0.  Box  209 

St.  Louis,  MO  63166 

Attn:  DRCPM-ASE  (J.  Keaton) 

Attn:  DRCPM-AS E-TM  (E.F.  Branhof) 

Attn:  DRCPM-AS E-TM  (MAJ  Schwend) 

Attn:  DRSAV-ASE-TM  (R.M.  Tyson) 

Attn:  DRSAV-EI  (CAPT  W.D.  Wolfinger) 

Attn:  DRSAV-EXH  (J.C.  Bulter) 

Army  Ballistic  Research  Laboratoreis 
Aberdeen  Proving  Ground,  MD  21005 
Attn:  DRXBR-VL  (R.G.  Bernier) 

Attn:  DRXBR-VL  (A.J.  Hoffman) 

Attn:  DRXBR-VL  (O.T.  Johnson) 

Attn:  DRXBR-VL  (R.  Mayerhofer) 

Attn:  DRXBR-VL  (D.M.  Mower) 

Attn:  DRXBR-VL  (W.S.  Tompson) 

Army  Materials  and  Mechanics  Research  Center 
Watertown,  MA  02172 

Attn:  DRXMR-R  (G.R.  Thomas) 

Attn:  DRXMR-RD  (R.W.  Lewis) 

Attn:  DRXMR-XC  (E.S.  Wright) 

Army  Materiel  Systems  Analysis  Activity 
Aberdeen  Proving  Ground,  MD  21005 
Attn:  DRXSY-J  (J.J.  McCarthy) 

Arnold  Research  Organization,  Inc. 

Arnold  Engineering  Development  Center 
Arnold  AFS,  TN  37389 
Attn:  Library 

Chief  of  Naval  Operations 
Washington,  D.C.  20350 
Attn:  OP-987  (Director) 

Defense  Documentation  Center 
Cameron  Station,  Bldg.  5 
Alexandria,  VA  22314 

Attn:  DDC-TRS-1 , 2 copies 

Department  of  Transportation  - FAA 
2100  Second  St.,  SW,  Rm.  1400C 
Washington,  D.C.  20591 
Attn:  ASD-520  (R.A.  Kirsch) 

Deputy  Chief  of  Staff  (AIR) 

Marine  Corps  Headquarters 
Washington,  D..  20380 

Attn:  AAW-61  (LT  COL  F.C.  Regan) 
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FAA/NAFEC 

Atlantic  City,  NJ  08405 

Attn:  ANA-r^O  (L.J.  Garodz) 

Attn:  ANA-64  (NAFEC  Library) 

Foreign  Technology  Division  (AFSC) 
Wright-Patterson  AFB,  OH  4543 
Attn:  FTD/NICD 

Marine  Corps  Development  Center 
Quantico,  VA  21134 

Attn:  D-042  (MAJ  W.  Waddell) 

Attn:  D-091  (LT  COL  J.  Givan) 

NASA  - Ames  Research  Center 
Mail  Stop  223-6 
Moffett  Field,  CA  94035 
Attn:  SC  (R.L.  Altman) 

Attn:  SC  (J.  Parker) 

NASA  - Ames  Research  Center 
Army  Air  Mobility  R&D  Laboratory 
Mail  Stop  207-5 
Moffett  Field,  CA  94035 

Attn:  SAVDL-AS  (V.L.J.  Di  Rito) 

Attn:  SAVDL-AS -X  (F.H.  Immen) 

NASA  - Johnson  Spacecraft  Center 
Housto,  TX  77058 
Attn:  JM6 

NASA  - Lewis  Research  Center 
21000  Brookpark  Rd. 

Mail  Stop  500-202 
Cleveland  OH  44135 

Attn:  Library  (D.  Morris) 

National  Bureau  of  Standards 
Building  225,  Rm  A62 
Washington,  D.C.  20234 
Attn:  I. A.  Benjamin 

Naval  Air  Development  Center 
Warminster,  PA  18974 

Attn:  Code  30C  (R.A.  Ritter) 

Attn:  Code  5422  (M.C.  Mitchell) 

Attn:  Code  5422  (C.E.  Mur row) 

Attn:  Code  5422  (D.G.  Tauras) 

Attn:  Code  5423  (B.L.  Cavallo) 

Naval  Air  Propulsion  Test  Center 
Trenton,  NJ  08628 

Attn:  ADI  (W.G.  Hawk) 

Attn:  PE3A  (J.  Mendrala) 
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Naval  Air  Systems  Command 
Washington,  D.C.  20361 

Attn:  AIR-0 3PAF  (CDR  R.C.  Gibson) 

Attn:  AIR-03PA4  (T.S.  Momiyama) 

Attn:  AIR-330B  (E.A.  Lichtman) 

Attn:  AIR-5 03W1  (E.A.  Thibault) 

Attn:  AIR-5204 

Attn:  AIR-5 204A  (D.  Atkinson) 

Attn:  AIR-5204J  (LT  COL  R.T.  Reisers) 

Attn:  AIR-53031  (R.O.  Lutz) 

Attn:  AIR-530313  (R.D.  Hume) 

Attn:  AIR-5323 

Attn:  AIR-53242  (C.F.  Magee) 

Attn:  AIR-5363 

Attn:  AIR-53632E  (C.D.  Johnson) 

Naval  Material  Command 
Washington,  D.C.  20360 

Attn:  MAT-0331  (H.G.  Moore) 

Naval  Postgraduate  School 
Monterey,  CA  93948 

Attn:  Code  57BP  (R.E.  Ball) 

Attn:  Code  57BT  (M.H.  Bank) 

Naval  Sea  Systems  Command 
Washington,  D.C.  20362 

Attn:  SEA-03511  (C.H.  Pohler) 

Naval  Ship  Engineering  Center 
Hyattsville,  MD  20782 
Attn:  Code  6105D 

Naval  Ship  R&D  Center 
Annapolis , MD  21402 

Attn:  Code  2851  (R.O.  Foernsler) 

Attn:  Code  2851  (J.R.  Lugar) 

David  W.  Taylor  Naval  Ship  R&D  Center 
Bethesda,  MD  20084 
Attn:  Code  1740.2  (F.J.  Fisch) 

Attn:  Code  1740.2  (O.F.  Hackett) 

Naval  Surface  Weapons  Center 
Dahlgren  Laboratory 
Dahlgren,  VA  22448 

Attn:  DF-52  (W.S.  Lenzi) 

Attn:  DG-10  (S.  Hock) 

Attn:  DG-10  (T.L.  Wasmund) 

Attn:  DG-104  (T.H.  McCants) 

Attn:  DK-23  (B.W.  Montrief) 

Attn:  DT-51  (J.F.  Horton) 
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Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
Silver  Spring,  MD  20910 
Attn:  WA-11  (E.F.  Kelton) 

Attn:  WU-41  (J.C.  Hetzler) 

Naval  Weapons  Center 
China  Lake,  CA  93555 

Attn:  Code  40  (M.M.  Rogers) 

Attn:  Code  40701  (M.H.  Keith) 

Attn:  Code  408  (W.T.  Burt) 

Attn:  Code  408  (H.  Drake) 

Attn:  Code  408  (C.  Padgett) 

Attn:  Code  4083  (G.  Moncsko) 

Naval  Weapons  Support  Center 
Crane,  IN  47522 
Attn:  Code  5041  (D.K.  Sanders) 

Pacific  Missile  Test  Center 
Point  Mugu,  CA  93042 

Attn:  Code  13^2  (J.R.  Bok) 

Attn:  Code  1332  (W.E.  Chandler) 

Attn:  Code  1332  (B.E.  Nofrey) 

Picatinny  Arsenal 
Dover,  NJ  07801 

Attn:  SARPA-AD-C  (S.K.  Einblnder) 

Rock  Island  Arsenal 
Rock  Island,  IL  61201 

Attn:  DRSAR-PPV  (D.K.  Kotecki) 

Attn:  DRSAR-RDG  (L.J.  Artioli) 

Attn:  DRSAR-SAS  (S.  Olsen) 

Attn:  SARRI-LR  (C.S.  Hicks) 

Attn:  SARRI-LW-W  (J.S.  Hansen) 

San  Antonio  Air  Logistics  Center 
Kelly  AFB,  TX  78241 
Attn:  ALC/MMSRE 

Warner  Robins  Air  Logistics  Center 
Robins  AFB,  GA  31098 

Attn:  WRALC/MMET  (LT  COL  G.G.  Dean) 

2750th  Air  Base  Wing 
Wright-Patterson  AFB,  OH  45433 
Attn:  ABW/SSL 
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Aeroqulp  Corp. 

Subsidiary  of  Libbey-Owens  Ford  Co. 

300  S.  East  Ave. 

Jackson,  MI  49203 
Attn:  E.R.  Steinert 

Attn:  R.  Rogers 

Air  University  Library 
Maxwell  AFB,  AL  36112 
Attn:  Librarian 

AVCO  Corp. 

Lycoming  Division 
550  So.  Main  St 
Stratford,  CT  06497 
Attn:  R.  Cuny 

Attn:  H.F.  Grady 

Battelle  Memorial  Institute 
505  King  Ave. 

Columbus,  OH  43201 
Attn:  J.H.  Brown,  Jr. 

Beech  Aircraft  Corp. 

9709  E.  Central  Ave. 

Whichita,  KS  67201 

Attn:  Engineering  Library  (T.R.  Hales) 

Bell  Helicopter  Co, 

A Textron  Co. 

P.0.  Box  482 
Fort  Worth,  TX  76101 
Attn:  J.R.  Johnson 

Attn:  J.F.  Jaggers 

Attn:  E.A.  Morris 

Boeing  Vertol  Company 
A Division  of  The  Boeing  Co. 

P.0.  Box  16858 
Philadelphia,  PA  19142 

Attn:  J.E.  Gonsalves  (M/S  P32-19) 

Cessna  Aircraft  Co. 

Wallace  Division 
P.0.  Box  1977 
Wichita,  KS  67201 
Attn:  B.B.  Overfield 

Attn:  Engineering  Library  (J.  Wilson) 
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Fairchild  Industries,  Inc. 

Fairchild  Republic  Co. 

Conklin  Street 

Farmingdale,  L.I.,  NY  11735 
Attn:  J.A.  Arrighi 

Attn:  Engineering  Library  (G.A.  Mauter) 

Falcon  Research  and  Development  Co. 

696  Fairmount  Ave. 

Baltimore,  MD  21204 
Attn:  W.J.  Douglass,  Jr. 

Falcon  Research  and  Development  Co. 

601  San  Pedro  NE,  Suite  205 
Albuquerque,  NM  87108 
Attn:  W.L.  Baker 

Fiber  Science,  Inc. 

7006  Sea  Cliff  Rd. 

McLean,  VA  22101 
Attn:  R.N.  Flath 


Firestone  Coated  Fabrics 
P.0.  Box  864 
Magnolia,  AR  71753 
Attn:  L.T.  Reddick 


Co, 


General  Dynamics  Corp. 

Convair  Division 
P.0.  Box  80877 
San  Diego,  CA  92138 

Attn:  Research  Library  (U.J.  Sweeney) 

Attn:  J.P.  Waszczak 

General  Dynamics  Corp. 

Fort  Worth  Division 
Grants  Lane,  P.0.  Box  748 
Fort  Worth,  TX  76101 

Attn:  P.R.  deTonnancour/G.W.  Bowen 


General  Electric  Co. 
Aircraft  Engine  Group 
1000  Western  Ave. 

West  Lynn,  MA  01905 
Attn:  E.L.  Richardson 

Attn:  J.M.  Wannemacher 
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General  Electric  Co. 

Aircraft  Engine  Group 
Evendale  Plant 
Cincinnati,  OH  45215 

Attn:  AEG  Technical  Information  Center  (J.J.  Brady) 

Goodyear  Aerospace  Corp. 

1210  Massillon  d. 

Akron,  OH  44315 

Attn:  J.E.  Wells  (D/959) 

Attn:  H.D.  Smith  (D/490G) 

Attn:  T.L.  Shubert  (D/910) 

Grumman  Aerospace  Corp. 

South  Oyster  Bay  Rd. 

Bethpage,  NY  11714 
Attn:  J.P.  Archey  (D/662) 

Attn:  R.W.  Harvey  (D/661) 

Attn:  H.L.  Henze  (D/471) 

Attn:  Technical  Information  Center  (J.  Davis) 

Hughes  Helicopters 
A Division  of  Summa  Corp. 

Centinela  & Teale  St. 

Culver  City,  CA  90230 

Attn:  R.E.  Rohtert  (15T288) 

Attn:  Library  (D.K.  Goss  2/T2124) 

ITT  Research  Institute 
10  West  35th  Street 
Chicago,  IL  60616 
Attn:  I.  Pincus 

Kamen  Aerospace  Corporation 
Old  Winsor  Rd. 

Bloomfield,  CT  06002 
Attn:  H.E.  Showalter 

Lockheed  - California  Co. 

A Division  of  Lockheed  Aircraft  Corp. 

Burbank,  CA  91503 

Attn:  Technological  Information  Center  (84-40,  U-35,  A-l) 

Lockheed  - California  Co. 

A Division  of  Lockheed  Aircraft  Corp. 

P.0.  Box  551 
Burbank,  CA  91520 
Attn:  C.W.  Cook 

Attn:  L.E.  Channel 
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Lockheed  - Georgia  Co. 

A Division  of  Lockheed  Aircraft  Corp. 

86  S.  Cobb  Drive 
Marietta,  GA  30063 
Attn:  C.K.  Bauer 

LTV  Aerospace  Corporation  / 

Vought  Systems  Division 
P.0.  Box  5907 
Dallas , TX  75222 
Attn:  G.  Gilder,  Jr 

Attn:  Unit  2-54244  (D.M.  Reedy) 

McDonnell  Aircraft  Co. 

McDonnell  Douglas  Corp. 

P.0.  Box  516 
St.  Louis,  MO  63166 
Attn:  R.D.  Detrich 

Attn:  R.A.  Eberhard 

Attn:  Library 

Attn:  M.  Meyers 

McDonnell  Douglas  Corp. 

3855  Lakewood  Blvd. 

Long  Beach,  CA  90846 

Attn:  Technical  Library,  (CL  290/36-84) 

New  Mexico  Institute  of  Mining  and  Technology 
Socorro,  NM  87801 
Attn:  TERA 

Northrop  Corp. 

Aircraft  Division 
3901  W.  Broadway 
Hawthorne,  CA  90250 

Attn:  Code  3680/35  (J.H.  Bach) 

Attn:  V.B.  Bertagna 

Attn:  Mgr.  Library  Services  (H.W.  Jones) 

Attn:  Code  3680/35  (W.  Hohlenhoff) 

Attn:  Code  3628/33  (J.R.  Oliver) 

Northrop  Corp. 

Ventura  Division 

1515  Rancho  Conejo  Blvd, 

Newbury  Park,  CA  91320 
Attn:  M.  Raine 
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Parker  Hannifin  Corp. 

18321  Jamboree  Rd. 

Irvine,  CA  92664 
Attn:  C.L.  Kimmel 

Attn:  J.E.  Lowes 

Protective  Materials  Co. 

York  St. 

Andover,  MA  01810 
Attn:  M.H.  Miller 

Rockwell  International  Corp. 

4300  E.  Fifth  Ave. 

P.0.  Box  1259 
Columbus,  OH  43216 

Attn:  Technical  Information  Center  (D.Z.  Cox) 

Rockwell  International  Corp. 

Los  Angeles  Aircraft  Division 
B-l  Division 
International  Airport 
Los  Angeles,  CA  90009 
Attn:  W.L.  Jackson 

Attn:  R.  Moonan  (AB-78) 

Attn:  W.H.  Hatton  (BB-18) 

Russell  Plastics  Tech. 

521  W.  Hoffman  Ave. 

Lindenhurst,  NY  11757 
Attn:  J.C.  Hebron 

Sikorsky  Aircraft 

A Division  of  United  Aircraft  Corp. 

Main  Street 
Stratford,  CT  06602 

Attn:  D.  Fansler/S.  Okarma 

Attn:  J.B.  Faulk 

Attn:  G.W.  Forbes 


Southwest  Research  Institute 
8500  Culebra  Rd. 

P.0.  Drawer  28510 
San  Antonio,  TX  78284 
Attn:  Bussuy-02 

Attn:  W.D.  Weatherford 
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Teledyne  CAE 
1330  Laskey  Rd. 

Toledo,  OH  43697 

Attn:  Librarian  (M.  Dowdell) 

Attn:  Librarian  (M.  Dowdell/W.Q.  Wagner) 

Attn:  Librarian  (M.  Dowdell/A. E.  Kirschmann) 

Teledyne  Ryan  Aeronautical 

2701  Harbor  Dr. 

San  Diego,  CA  92112 

Attn:  Technical  Information  Services  (W.E,  Ebner) 

Attn:  P.  Kleyn 

Attn:  N.S.  Sakamoto 

The  Boeing  Co. 

Wichita  Division 
3801  S.  Oliver  St. 

Wichita,  KS  67210 

Attn:  H.E.  Comer  (M/S  K21-57) 

Attn:  D.Y.  Sink  (M/S  K16-14) 

Attn:  Library 

The  Boeing  Co. 

Aerospace  Group 
P.0.  Box  3999 
Seattle,  WA  98124 

Attn:  J.G.  Avery  (M/S  41-37) 

Attn:  R.G.  Blaisdell  (M/S  8C-42) 

The  Johns  Hopkins  University 
Applied  Physics  Laboratory 
Johns  Hopkins  Road, 

Laurel,  MD  20810 
Attn:  A.R.  Eaton 

Attn:  C.F.  Meyer 

Attn:  B.W.  Woodford 

Uniroyal , Inc . 

Mishawaka  Plant 
407  N.  Main  Street 
Mishawaka,  IN  46544 
Attn:  J.D.  Galloway 

Uniroyal , Inc . 

Government  Affairs 
1700  K.  St. , NW 
Washington,  DC  20006 
Attn:  D.  Gillett 

Williams  Research  Corp. 

2280  W.  Maple  Rd. 

Walled  Lake,  MI  48088 
Attn:  Library 


£ S 8ai 

.'HJi 

: « «»  ^ S 

fi'fr 

:T  1 1 1 
: c u "s  i 
' ^ *;  C o 

«o 

0 2 £ °r  w 

: £ £ £ 15  ! 

0-  a< 

! ! c 

1 # c sr'S 

' _ § a 

l SQ  i S ! 

a.  : i>  <— 

i'S**  . 

; * ^ H a-  • 

■■  .£  o U.  . 

i s < 

; t/i  SP  o 

: o *■  i 

S S O X 

; 5 Q o O 1 

* *2  .!£}  Q> 

s ^ > a 
(jSOm 


i-  LU  £ 
O'  — — 
— LU  3 

a 

E«  2 

£y  ~ 

,*j22 

c 

>,  c » 

4->  tO  CD 

3?  & 

•2  o 

.ass 

3 a. 
5 cl  57 


2§F 

o > ^ 


O w)  n 

- Sw 

- £ * O 

1)1  I* 

t?  -a  c 
>.  o « 
1/1  E £ 

7 ! e/i 

If  » 

'C  U c 

I 2 « ® 
i u t:  n 
i u.  C 

*3  J2  -tS 
c ** 

;J!<  o 

! I S* 
; ;=  &£ 

! is  2 .2 

,25£ 

! c Jj  & 

*3  *£3 

! 2 g £ 

i M f s 

[ *«-*  CQ 
. -C  ~ 
v-  —j 

L O "3  JJ 

: « % S 
L-S  S-o 
: 

}•§  I & 

■ E x:  O 
’ « ~ " 

^ ti  > 

>>  jo  <d 
•a  *3  T3 


o 


J3,  s>  O CL 1 

“ te  o ■ 

«■:  o ^ J ( 

:*2  «- ^ j 

i'i<.r 
* c § 
o e •p  i 


2 o 

O w 


51^3 

sSug  5 

&J  o 
,x> 

£ "O  H 


# g sl  'O 

■ .ig^ 

?i°s£ 

■ 5»  oi 

^ 5f  ,_•  J 

^ H ft, 

j£<  < 

• .£  o U. 

i t3  --  < 


. 5“  - I 

; :=-o  § o 

; | < I » 

CsDffl 


cL  ^ 

u « 
n 

sO  ^ •— 

o'  d ° 

r~  uj  « 

— E 5 

St« a 
£<  - 

iy 

uZ  a 
to  D c 

^§! 
i’s  s. 
•s  ^ - 

2 o 
- a 

3 a ^ 

4-  ^ </> 

E § p 

o > ^ 


0 ifi  « 

w " , 

, g 5 © 

4/  <u  <2 

y>  "2  c 
X o u> 

" E £ 

1 og1 
c;  *5  "» 

•c  3 c 

2 «5 
S c w 

”IS 


> c .t: 
COS 
33  O.X) 

• — rt  sj 
O u £ 

2 5o 
g a»  a 

(A  pC 

v a> 

•£  >.« 
v-  ^ 3 

o TJ  « 

- K 2 
^ o u 
o °«  o 

F in" 

~ « "O 

■a  K 21 

£ o 

« *-  u 
= w > 
•a  ■»  -a 


O 


WWW«W8WW»»aB^^ 


(Contd.  on  Card  2) 


o <*-• 

g e 

•a  <o 
53  £ § 

CJ  u- 

voE 

: gS 

: i * 

• <D  P.4 

§?£l 

c~  -a 

fti  <1>  _* 

•I  £ * 


>§  2<~ 
a o 

2 fc  | >. 
2 & * 
O 3 S ~ 

x>  ft.  « x> 

2 5 g 2 
to  .|  g 

c VJ  St'S 

O >s  " 3 

•g  o *3  > 

3 >v*2  <u 

2 xi 

o ? E“ 

u.  -w 

fatO  MC 
o >,.S  ° 

u 

s 2 3 a 

< •»  v>  ca 

« S £ ts 
2=5  «- 
o 2 x « 

u-i;Pa 

3 ‘s 


5/5 

8 * 

<u 

J8  6 


6 c . 

y 9>  T3 
C > « 
•—  l-  T3 
*j  O c 
O . g 
fisc 
S.2E 
8 3 8 

^ J 2 


>>  x £ 

jS  8 I 
2 » 1 
c«  — « g< 
J-  V O 

ja  « 

3^  <u 
> O T3 


o 


■«-»  M*  c/5  'TJ 

&g<£  § 

, c c G +i 

•S'3ug 
£ g 

^ 4-h  .ia 

o E o 

2 ~ g £ 

5 >,.2  8 

a «?  *3  J3 

TO  TO  4—i 

« ^ JH  t 

> c <1)  Ji  3 

’ e a 3 fcu 
•o  «»  • 

S?  « C b 

8-hu? 

S’  ° £ S 
S >,  2 ® 

c 3 a 60 
« 2 £?« 

* ’’33  O G 

ft  2 O *>  T3 
§J  a>  e ^ w 
.S3  c *3  T3 
»3  +-  o c 

.§  s § c e 

'?£  s-S  | 

2 - o g © 

M C ^ x 2 

C o +■>  T* 

3 .22  >.x  ■£ 
SSSoS 

U t/5  »PH  pH  C 
n tfl  C 

£i  2 y §■ 

L 3 t>  >h 

^ 3*2  ^ 


Jii)ilWW^i!M»W 


